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TRAVELING WAVE FOCUSING FOR PLASMA CONTAINMENT 


C. K. Birdsall and Allan J. Lichtenberg 
Department of Electrical Engineering, University of California, Berkeley, California 
(Received June 29, 1959) 


It is a well known fact that particles moving 
through a system of periodic lenses are focused.’ 
Atraveling wave, when viewed in a system mov- 
ing with the wave, appears as a system of elec- 
tric and magnetic lenses. If a relatively station- 


ary plasma is viewed from this moving coordinate 


system, the charged particles, both electrons 
and ions, move through the lens system at the 
phase velocity of the wave. Both electrons and 
ions are focused in such a system with the solu- 
tions given in the paraxial approximation by the 
Mathieu equation. One may equally well analyze 
the system in the laboratory coordinates with 
moving lenses and obtain the same result. For 
magnetic fields alone the equation of motion can 
easily be shown to be 


dy lw? 
doh? * 8 Gr (1+ cos2w#)r =0, 


where ry is the radius, w the applied frequency, 
and w, the cyclotron frequency corresponding to 
the peak magnetic field. There are a series of 
focus and nonfocus bands depending on the values 
of the applied magnetic field and the applied fre- 
quency. Also, particles traveling synchronously 
with the wave see a constant field and should not 
be focused. 

Both of these methods of losing particles seem 
tobe absent if a periodic structure with spatial 
harmonics is used. The waves associated with 
the harmonics have different power levels and 
for particles moving with a finite velocity also 
different apparent frequencies. The apparent 
frequency of variation of the fields as seen by the 





wave is 


“sop =W- Bl 
where £,,, is the propagation constant of the mth 
harmonic and u the axial velocity of the particle 
in question. 

In the experiments described a low-density 
plasma was used. If the plasma density is high 
the fields will be excluded from the plasma, but 
focusing of the same nature should take place in 
the high-gradient fields near the edge of the 
plasma. In this case the energy density of the 
fields must be greater than the energy density of 
the plasma for containment. 

Two experiments were performed with plasma 
columns extending through slow wave structures 
3 cm and 23 cm long, respectively. In both cases 
the diameter of the plasma column was 0.75 cm. 
Focusing was tried in the frequency range of 3 to 
25 Mc/sec and at 1000 Mc/sec. In the low- 
frequency range a column of tens of milliamperes 
discharge could be focused by some tens of watts 
of power at all frequencies. The amount of con- 
traction of the plasma appeared proportional to 
the applied power over the entire frequency 
range. The focusing occurred at power levels at 
which the oscillation amplitudes of the electrons 
in a magnetic field as calculated from approxi- 
mate solutions to the Mathieu equation lay within 
the tube. There was also evidence of plasma ex- 
citation and heating in the increased brightness 
of the plasma with the applied rf and the extent 
that the plasma was extended in the longitudinal 
direction. The excitation of the plasma would be 
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caused by the presence of the electric fields. 

The focusing could still be observed at 100 ma 
arc current but did not appear as marked. In 
these experiments it is believed that the electrons 
are being focused by the traveling fields and the 
ions are contained in the potential well thus 
formed. With large structures and higher power 
it should be possible to focus both electrons and 
ions. At 1000 Mc/sec the ionization effects of the 
high-frequency effects were very strong and the 
plasma could be maintained without a dc dis- 
charge. No focusing was observed, but this may 
have been masked by the intense excitation. 

From the experiment to date one cannot say for 
certain whether the focusing was due primarily to 
electric fields or primarily to magnetic fields. 
The slow wave structures were designed to con- 
centrate the magnetic fields in the vicinity of the 
plasma but the observed excitation indicates that 
electric fields were also present. Other observa- 
tions also seemed to indicate the presence of 
rather strong electric fields. 

The circularly symmetric structures used con- 


fined the plasma in the transverse direction but 
not at the ends. A re-entrant plasma column 
(e.g., toroid or figure 8) may be focused using a 
re-entrant circuit excited to have only a traveling 
wave. Excitation may be achieved by the focusing 
field itself. In this and most other forms the cir. 
cuit power may be recirculated supplying only the 
losses in the structure and the plasma. 

The reader should note the analogy between the 
traveling waves discussed by Wuerker, Shelton, 
Langmuir and Goldenberg.’ 

The authors would like to acknowledge the in- 
teresting and helpful discussions they had with 
G. Kino (of Stanford University) on the general 
subject of plasma confinement. Dr. Kino and Dr, 
Chodorow were working on a quadrupole rf sys- 
tem at the time we conceived of our confinement 
scheme. 





1A. M. Clogston and H. Heffner, J. Appl. Phys. 25, 
436 (1954). ™ 
*Wuerker et al., J. Appl. Phys. 30, 442 (1959); 

Wuerker et al., J. Appl. Phys. 30, 441 (1959). 





We wish to report the first observation of nu- 
clear resonance in a ferromagnet. A remarkably 
strong resonance was detected at room tempera- 
ture in a sample of finely divided face-centered- 
cubic cobalt (Co**) metal.’ The resonance fre- 
quency was measured to be 213.1 Mc/sec in zero 
applied field. The observed frequency indicates 
a magnetic field of 213 400 oersteds at the site of 
the cobalt nuclei. This very sizable magnetic 
field arises from both core and conduction elec- 
tron contributions as has been discussed by 
Marshall.? 

The above results are in good agreement with 
the hyperfine field deduced from specific heat 
measurements on bulk hexagonal cobalt and 
cobalt alloys between 0.3 and 1.0°K.*»* The spe- 
cific heat work indicates a nuclear field of 220000 
oersteds in hexagonal cobalt. This work also 
suggests that the hyperfine field in the cubic 
phase is about 10000 oersteds lower than in the 
hexagonal phase. Only a minor temperature var- 
iation in the nuclear field is expected between 
0°K and room temperature. Thus the resonance 
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measurements and the specific heat measure- 
ments are in good agreement. 

A tracing of the resonance line is shown in 
Fig. 1. The half width at half maximum absorp- 
tion is 275 kc/sec. We believe that this line 
width is associated with a static broadening 
mechanism rather than with relaxation of the 
cobalt spins. In some respects the line appears 
to be inhomogeneously broadened. The observed 
line shape is that of an absorption curve rather 
than the sum of absorption and dispersion deriva- 
tives as is usually observed in metallic reso- 
nances. The very similar behavior of the F- 
center resonance in alkali halides has been esta- 
blished as an adiabatic rapid passage effect.* 
There are differences in detail however in the 
dependence of the resonance signal on modulation 
amplitude and frequency. The present signal is 
strictly proportional to the modulation amplitude 
over the range of frequencies investigated. At 
the higher modulation frequencies a decrease it 
amplitude and a phase shift develop. The ampli- 
tude and phase of the resonance signal as a func: 
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FIG. 1. Recorder tracing of nuclear resonance 

signal in finely divided fec Co®® metal. The oscillator 

%, is frequency modulated at 1.2 kc/sec with an ampli- 
tude of 200 ke/sec. 
tion of modulation frequency are shown in Fig. 2. 
The theoretical curves are for exponential re- 
laxation with a relaxation time of 0.13 millisec- 
onds. This time is somewhat shorter than the 
relaxation time expected from interaction with 
conduction electrons.® Still, the interaction of 
the nuclear spins with spin waves is expected to 
have an even smaller effect.” There exists the 
possibility that the observed relaxation time is 

P associated with spin-spin interaction. A strong 
coupling of the cobalt nuclear spins through spin 
waves can be expected.® 

FP We would next like to comment on the remark- 
able strength of the resonance. Ordinarily only 
avery small absorption might be expected from 

a samples with line widths as large and radio- 
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frequency penetration as low as in cobalt metal. 
However, built into the ferromagnetic sample is 
a mechanism which increases the strength of the 
nuclear resonance absorption by four orders of 
magnitude. It is this mechanism which is re- 
sponsible for the tremendously enhanced signal 
which we have observed. It can be seen that the 
resonance is excited predominantly by the hyper- 
fine coupling itself rather than by the external rf 
field and this leads to the enhancement. Since 
the magnetization relaxes rapidly, it can be pre- 
sumed to follow the rf field. A transverse field 
H, should produce an internal field H,/(1+ Nx) 
where N is the demagnetizing factor. We can 
then expect the magnetization to rotate through 
an angle H,/(1+Nx)H A Where H, is the anisotropy 
field. The rotation of the magnetization at the 
radio frequency should produce a transverse rf 
field of magnitude [H ,)/(1+ NX)H ,]H, where Hy 
is the hyperfine field at the cobalt nuclei. We 
estimate the factor H,/(1+NX)H, to be of the 
order of 107. Then the resonance intensity which 
is proportional to this factor squared should be 
enhanced by some four orders of magnitude as 
observed. One would expect that this enhance- 
ment factor would decrease in a static magnetic 
field as shown in Fig. 3. For comparison we 
have also measured the field dependence of the 
transverse susceptibility of our sample. One 
would expect that at low fields both domain wall 
displacement and domain rotation would contri- 
bute to the transverse magnetization. At high 





tT T T i 


a a a 


1 





© NUCLEAR RESONANCE 
ABSORPTION 

m= SQUARE OF TRANSVERSE 
SUSCEPTIBILITY 


ABSORPTION 
2 
. vv . S? TT 


— =a wee! 


° 


T 
L 








0.01 mame net) ! Ll 
100 1000 


FIELD (Oersteds) 





FIG. 3. Comparison between the strength of the 
absorption signal and the square of the transverse 
magnetic susceptibility as a function of static applied 
field. 
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fields only domain rotation should make an ap- 
preciable contribution. One can consider that 
the resonance signal measures the domain- rota- 
tion contribution to the magnetization. From the 
fact that the resonance signal drops off much 
more rapidly than does the square of the suscept- 
ibility we must conclude that domain wall pro- 
cesses play an important role out to the highest 
fields in which the resonance could be detected. 
The resonance signal has been observed in 
static magnetic fields up to 5400 oersteds. Above 
5400 oersteds the signal intensity drops below 
the noise level. Over this range of applied fields 
the resonance frequency holds constant to within 
100 kc/sec. This is somewhat surprising since 
one might expect a field of 5400 oersteds to pro- 
duce a frequency shift of some 5.4 Mc/sec. How- 
ever the cobalt particles because of their high 
demagnetizing fields are still relatively perme- 
able up to the highest fields at which the reso- 
nance could be observed. A plot of the sample 
magnetization obtained by integrating the apparent 
longitudinal susceptibility is shown in Fig. 4. At 
5400 oersteds the magnetization is still only 70% 
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FIG. 4. Magnetization vs applied magnetic field for 
finely divided fec cobalt as obtained from an integra- 
tion of the longitudinal ac susceptibility. 











saturated. We must conclude that variations in 
internal field are kept below a few hundred oer- 
steds over the range of fields investigated. 

The measurements reported in this Letter 
were made by placing the sample at the shorted 
end of a quarter-wavelength coaxial cavity. The 
cavity forms the resonant element of a 464A 
lighthouse -tube oscillator. Our spectrometer is 
very similar to a unit developed independently by 
Kojima et al.° The spectrometer is frequency- 
modulated by vibrating a small capacity across 
the high voltage end of the cavity. Modulation 
frequencies from a few cycles up to 10 kc/sec 
and modulation amplitudes up to 3 Mc/sec are 
possible with this design. Further studies both 
of the static and dynamic interaction between nu- 
clear and electron spins in a ferromagnet are in 
progress. 

We wish to acknowledge very helpful conversa- 
tions with Dr. W. Marshall, Dr. P.-G. de Gennes, 
and Professor C. Kittel. Our very deep gratitude 
is due Dr. Nicholas Kurti who aroused our in- 


terest in this problem during his stay in Berkeley. 


We are indebted to Mr. Harry Robinson for ad- 
vice on the construction of the spectrometer. 
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'We are indebted to Dr. J. Crangle for making this 
sample available to us. 
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DIRECT OBSERVATION OF PHONONS DURING TUNNELING IN NARROW JUNCTION DIODES 





N. Holonyak, Jr., and I. A. Lesk 
Advanced Semiconductor Laboratory, General Electric Company, Syracuse, New York 


and 


R. N. Hall, J. J. Tiemann, and H. Ehrenreich 
General Electric Research Laboratory, Schenectady, New York 
(Received July 14, 1959) 


Esaki has shown that narrow p-n junctions be- 
tween degenerate n- and p-type regions of a semi- 
conductor can give rise to a pronounced current 
maximum in the direction of forward bias.’ This 
component of the forward current is attributed to 
the tunneling of electrons from states lying near 
the bottom of the conduction band to states near 
the top of the valence band. Keldysh and Kane?»® 
have shown that in a semiconductor such as Ge or 
Si, where the conduction and valence band edges 
lie at different points in k space, the tunneling 
process may require the participation of one of 
the phonons having a wave number corresponding 
to the separation of the two minima. These pho- 
nons are the ones involved in optical transitions 
inthese semiconductors. This paper reports 
measurements on tunnel diodes made from ger- 
manium and silicon which show distinct structure 
at voltages corresponding to the onset of these 
phonon-assisted processes, and which permit a 
precise measurement of the energies of the pho- 
nons involved. In addition, the experiments pro- 
vide further insight into the relative probabilities 
of direct and indirect tunneling processes. 

The phenomenon which we report was first 
observed by Holonyak and Lesk in heavily doped 
silicon junctions, cooled to 4.2°K. Distinct dis- 
continuities in slope (cusps) of the current-voltage 
characteristic were noted for small forward and 
reverse voltages. The lack of sharpness of these 
cusps is attributed largely to thermal broadening. 
The current-voltage characteristic of one of 
these diodes is shown in Fig. 1. This junction 
was constructed with such a small donor concen- 
tration that the negative-resistance region is 
nearly absent. The cusps were found to occur in 
all of the silicon junctions which exhibited a meas- 
uable tunnel current. They appeared at the same 
voltages, for junctions made using P, As, or Sb 
asthe donor impurity, and did not vary as the 
donor concentration was changed in the range 
tetween 2x 10° and 10” cm~*. They did not de- 
pend upon crystal orientation, and were clearly 
observed using polycrystalline material in which 
the junction covered several different grains. 


~ ARBITRARY SCALE 


CURRENT 














FIG. 1. Recorder tracing of current-voltage char- 
acteristic of a Si junction at 4.2°K. The cusps at the 
lowest energies are identified with transverse acousti- 
cal (TA) and probably optical (O) phonons. The third 
cusp, which corresponds to a phonon or phonons ob- 
served in optical data, remains unidentified. 


Similar cusps were observed at 4.2°K in tunnel 
diodes made from germanium when Sb was used 
as the donor impurity, as shown in Fig. 2. How- 
ever, they appear much less prominently in diodes 
made with P or As as the donor impurity, using 
a range of donor concentrations (0.8 x10" to 
6 x10'® cm™’) which included that used for the Sb- 
doped diodes. 

We interpret the cusps as evidence that lattice 
phonons are required in order to conserve momen- 
tum during the tunneling process. The cusps occur 
in pairs, at equal forward and reverse voltages, 


and their locations agree within experimental error 


with the phonon energies which are deduced from 
optical measurements of Si and Ge.*»> In Ge, we 
find these energies to be (7.9+0.1 and 27.5+0.5) 
x10~* ev, respectively, and it is clear that they 
correspond to the transverse and longitudinal 
acoustic phonon energies. At present, we are 
unable to resolve clearly any structure which can 
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FIG. 2. Recorder tracing of current-voltage char- 
acteristic of a Ge junction at 4.2°K. Arrows indicate 
the locations of the cusps which are identified with the 
transverse and longitudinal acoustic (TA and LA) pho- 
nons. 


be attributed to optical phonons. In Si, three 
cusps are observed at (18.0+0.5, 5741, and 
120+5)x10™ ev. The first corresponds to the 
transverse acoustic phonon energy, and the sec- 
ond appears to be due to one of the optical pho- 
nons. The third is unidentified, although it is 
also observed in the optical data.*»® 

We examined tunnel diodes made from GaSb 
and found no evidence of any cusps. This behav- 
ior is understandable inasmuch as both band 
edges are located at k =0 in this material,°® so 
that tunneling can proceed without assistance 
from phonons. 

In germanium junctions made using P or As as 
the donor impurity, the cusps are difficult to re- 
solve, indicating that the tunnel current in these 
junctions results largely from transitions which 
do not involve lattice phonons. This dependence 
upon the specific donor impurity is difficult to 
understand, and will be the subject of further 
study. It is conceivable that with large concen- 





trations of these impurities, the density of states 
in the conduction band is modified in such a way 
that an appreciable number of conduction band 
states are found at k =0 below the Fermi level, 
and that the degree of this modification depends 
specifically upon the impurity involved. In pure 
germanium, the k=0 minimum is only 0.15 ey 
above the (111) minimum,’ and therefore this 
heavily doped material might have its band edge 
at (000) modified sufficiently to permit an appre- 
ciable rate of direct tunneling transitions. This 
tailing of the band edge is similar to that calcu- 
lated by Parmenter® who considered the effect of 
moderately large concentrations of hydrogenic 
impurities on the band structure in the vicinity 
of a simple edge. Additional evidence for tailing 
of the band edge is provided by studies of the 
shape of the negative-resistance region in tunnel 
diodes. It is found that this region always extends 
to higher voltages than is expected from the 
known impurity concentration and the density of 
states of the pure semiconductor. 

The authors would like to thank R. L. Watters 
for assistance in obtaining the data; G. K. Wessel 
and J. Young for making their liquid helium facili- 
ties available; and S. F. Bevacqua for assistance 
in developing suitable materials for fabricating 
silicon tunnel diodes. 
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MECHANISM FOR ANOMALOUS MAGNETOCRYSTALLINE ANISOTROPY PEAKS 





IN FERROMAGNETIC CRYSTALS* 


C. Kittel 
Miller Institute for Basic Research in Science, University of California, Berkeley, California 
(Received July 20, 1959) 


Dillon’ reported the occurrencé of strong sharp 
peaks in the magnetocrystalline anisotropy con- 
tribution to the ferromagnetic resonance field in 
yttrium iron garnet (YIG) at liquid helium tem- 
peratures. His results have been confirmed by 
Dillon and Nielsen,? who have attributed the anom- 
alous peaks to rare earth ion impurities, and 
striking results for Tb-doped YIG were presented. 
At 1.5°K at 23 kMc/sec in YIG containing 0.1% Tb 
three or more anomalous peaks are observed on 
rotation in the (110) plane, the height (always 
positive) of the peaks being of the order of 3 to 5 
koe and the angular width at half-maximum of the 
order of 3 deg. The sharpness of the peaks defy 
reasonable analysis in spherical harmonics, and 
no previous theoretical attempt to explain the 
peaks has broached the problem, although Wolf® 
has recognized that the exchange coupling of Fe** 
and rare earth ions must be involved, and Dillon 
and Nielsen? state that Clogston, Walker, and 
Dillon are attacking the problem. 

The underlying mechanism appears to be very 
simple, and we give here an elementary model 
which makes the physics clear. When the crystal 
field splitting parameters become known for the 
relevant rare earth impurities in YIG or in the 
gallate, it will be easy to make detailed predic- 
tions. Consider as in Fig. 1 the energy levels of 
a paramagnetic ion in a uniaxial crystal field and 
an effective exchange field arising from the ferric 
lattice magnetization at an angle 6 with the crys- 
tal field axis. For the highly anisotropic g’s 
characterizing most rare earth levels in crystals 
it is legitimate to treat the exchange interaction 
as an effective magnetic field.* To obtain an 
anomalous anisotropy it is necessary that the 
ground states cross over® at one or more angles 
§. At absolute zero the ion is in the state m, 
when the magnetization lies between cos~0 and 
cos*@,, and in the state m, between cos~*é, and 
cos1. The anisotropy field H, in a resonance 
experiment is essentially defined by 


H =- #E/ae 
a”s E/3@, (1) 
where M, is the total magnetization and E refers 


to unit volume. At absolute zero H,, will have a 
delta -function type maximum at each crossover 


ENERGY 











Eg---— 
| ” 
= L cos @ 
0 6 t 
FIG. 1. Schematic crystal field level diagram in com- 


bined exchange and uniaxial crystal fields, where @ is 
the angle between the crystal field axis and the direc- 
tion of the effective exchange field arising from the 
ferric lattices. 


6. On either side of 6, there will be a weak 
anisotropy field of the order mH ext cosé, where 
the effective quantum number or g-value m is 
m, OF m,, according to the direction, and f is 
the ratio of rare earth ions to unbalanced ferric 
ions; for f ~10~* and® H,, ~ 10° this background 
anisotropy of ~100 oersteds from the rare earth 
impurities is to be added to the sources of ani- 
sotropy in pure YIG. 

The delta-function anomaly spreads out at a 
finite temperature. We consider deviations Aé@ 
from 6, such that (Am)uH,,40« kT. Within this 
angular range we see from the equilibrium popu- 
lation difference that 


(AE) =~ (Am)*(uH , 46)*/kT, (2) 


within a factor of the order of unity. Thus the 
anomalous anisotropy field is given by 
H ~ (Am)*(uH *)f /kT = 10*/T oersteds, (3) 


where T is in deg K. The anomaly has an angular 
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laxation processes are not essential. Exchange 
interactions among the rare earth ions and in- 
directly via the ferric ions will also limit the 
temperature to which our results apply. These 
systems may have important applications to mag. | 
netic cooling by turning the external field toa 
crossover point. 

I am indebted to P.-G. de Gennes for helpful 
discussion and to W. P. Wolf for a preprint in 
advance of publication. 





width 
Aé wkT/uH Am =0.037T radian. (4) 










These estimates are reasonable agreement with 
the observations. The order of magnitude of 
height and width are correct, and the directions 
of the temperature dependences are correct. It 
should be noted that our mechanism predicts the 
anomaly should always be an increase in reso- 
nance field, as seems to be observed. Different 
ion sites will give crossovers at different angles. 
The model makes one important nontrivial 
assumption —that at the crossover point the rare 
earth ions relax rapidly so that even at micro- 
wave frequencies the ions stay in thermal equilib- 
rium. Otherwise we would lose the anomaly. A 
phase lag will lead to a line broadening at the 
anomaly. The lines are indeed considerably 
broadened at the anomaly, but further experi- 
ments are required to determine the actual mag- 
nitude of this possible contribution to the width. 
If, however, the levels do not actually cross over, 
but are split by perturbations at the nominal cross- 
over, then Eqs. (3) and (4) are valid only for 
T>T,=4/k, where A is the splitting. Now re- 
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In previous experiments originally undeformed, thin aluminum foils, containing quenched-in dis- ¥ 

electron-transparent aluminum foils have been location loops, while being strained in the elec- 2 
mechanically strained while under observation tron microscope. It was found that in them the 3 
in the electron microscope.’ Among the inter- typical smooth dislocation lines arose which 
resting results, published elsewhere,” was the often moved quite closely past loops without 
observation that the dislocation lines introduced obviously affecting them. In other cases loops 
during the straining are mostly smooth and move moved out of the specimens, either spontaneously, 
without drastically changing their shape. From or under the influence of an approaching glide 
other work it is known that rapidly quenched dislocation; and occasionally small loops were 
aluminum contains a high concentration of pris- observed to emit glide dislocations when a dis- 
matic dislocation loops, due to the condensation location moved past them at a close distance, # 
of thermal vacancies,* as had been predicted effect which was predicted in an earlier paper.’ F 
previously.*; ° In an effort to determine what influence, if a 


Recently a study was made of the behavior of 
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any, specimen size has on these phenomena, 
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quenched aluminum was deformed approximately aggregates would be invisible in the electron 
2%, before thinning it into electron-transparent microscope, but dislocations would be held up 
foils. The surprising result was that, instead of when and where their cores meet an aggregate. 

the smooth dislocation lines, intricate tangles of The vacancies would then begin to diffuse into 
kinked and branching dislocations with inter - and along the cores to form or add to an already 
spersed loops were found (Fig. 1), whereby the existing atmosphere of vacancies around the dis- 
loop density seemed to be somewhat smaller location axes. Once a sufficient vacancy concen- 
than in the as-quenched specimens. tration has thus been reached condensed vacancy 

Finally, furnace cooled aluminum was lightly disks should be nucleated, aided by the stress 
strained, and then thinned for observation. Even field of the dislocations. If the new prismatic 
this material showed similar dislocation struc - dislocation loops are formed so as to partly fuse 
tures, but by far fewer and on a coarser scale, with the original dislocations they form kinks. 
(Fig. 2). Such kinks either pin the dislocation entirely or 

How may we interpret these findings? It is are dragged along with them, or they may be 
suggested that quenched aluminum —and to a “shaken off,” that is they may be left behind in 
lesser extent furnace cooled aluminum — contains the form of closed, irregularly shaped loops. 
ib- vacancy aggregates of between a few and a few If, on the other hand, the condensation disks 
hundred vacancies each, an idea which agrees are formed separate from the dislocations nu- 

. with findings of Federighi and Panseri.’ Such cleating them, then they will be simply left be- 
Phys hind as the dislocation moves on. In a similar 
fashion small or unfavorably oriented dislocation 
loops in the path of a moving dislocation would 
be dissolved in favor of generating kinks or 
larger or more favorably oriented loops. Further, 
kinks and loops themselves may act as disloca- 
tion sources, particularly, as we saw, in the 
vicinity of other dislocations. The result of all 
these processes might be described as a “mush- 
rooming” of any dislocation starting to move, 
yielding the observed dislocation structures of 
Figs. 1 and 2. 

In the very thin foils, as may be investigated 
by transmission microscopy, any vacancies in 
atmospheres around dislocations, however, will 

FIG. 1. Transmission electron micrograph of diffuse out of the specimen surfaces, so that the 
quenched and subsequently deformed aluminum. 40 000:1 supersaturation necessary for the nucleation of 
collapsed vacancy disks is never attained. The 
same must hold true for comparably thin surface 
layers in thick specimens. 

Vacancies generated through slip should give 
rise to the same effects. Thus, moving disloca- 
tions would generate vacancies, then allow them 
to precipitate into loops or kinks which partly 
block the progress of the original dislocations, 
but also act as dislocation sources, greatly in- 
creasing the dislocation density. Work-hardening 
in metals may be partly due to this phenomenon, 
particularly since the described dislocation 
“tangles” are not only found in aluminum but 
‘ have also been observed when repeating the ex- 
2, at _—— a periment with nickel. At very low temperatures 

FIG. 2. Transmission electron micrograph of fur- of deformation, interstitials would assume the 
nace-cooled and subsequently deformed aluminum. role played by vacancies at higher temperatures, 
40000:1. except that they condense only in the form of 
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disks but not of three-dimensional aggregates. 

In many alloys, on the other hand, vacancies 
are bound to the substitutional atoms, and condi- 
tions would be completely different from those 
discussed here, but the profuse multiplication of 
dislocations in crystals originally almost free of 
dislocations might in some cases be due to the 
phenomena discussed. 





This research was supported by the U. S. Atomic 
Energy Commission. 
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NATURE OF RADIATION DAMAGE IN NICKEL* 


H. G. F. Wilsdorf 
The Franklin Institute Laboratories, Philadelphia, Pennsylvania 
(Received June 18, 1959) 


The problem of radiation damage by neutrons 
has been reviewed recently by various authors.’»? 
The evidence for the damage in metals obtained 
by resistivity measurements, x-ray techniques, 
tests of mechanical properties, and others can- 
not be interpreted unambiguously. In order to 
attempt a better understanding of the processes 
involved, a more direct experimental approach 
was sought by investigating irradiated metal 
foils by transmission electron microscopy. This 
method has been used for directly studying dis- 
locations, stacking faults,°~* and dislocation 
loops produced by the condensation of vacancies® 
Also, it is a most sensitive tool to detect lattice 
misorientations which may be as small as 1 
minute of arc. The disadvantage of the method 
is that, at present, it does not permit tests be- 
low room temperature. However, since many 
metals show a considerable increase in yield 
strength after neutron exposure at pile tempera- 
ture’~® the described experiments appeared to be 
suitable for investigating the nature of radiation 
damage that can be produced in this temperature 
range. The tests were made with nickel speci- 
mens of 99.999% purity which had been irradiated 
as foils 1000 A to 2000 A thick as well as with 
thicker specimens. The latter had to be thinned 
after radiation treatment for their examination. 
The transparent areas were obtained by electro- 
lytic polishing from rolled strip. 

Stimulated by the theoretical work on displace- 
ment spikes,’® and by the theory that possibly 
very small dislocation loops might be produced 
by a knock-on," foils in the thickness range of 
1000 A were exposed to neutrons. All specimens 
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were sealed in evacuated (10°° mm Hg) quartz 
tubes and irradiated in the Brookhaven reactor.” 
The same areas were examined at magnifications 
between 40000x and 65 000x before and after 
irradiation with an integrated flux of 10° nvt or 
10** nvt. Even the most careful scrutiny of the 
many specimens investigated failed to reveal any 
consistent structural changes in the foils. At 
this point it seems to be safe to conclude that 
the neutron bombardment of thin nickel foils at 
pile temperature does not produce dislocations 
or stacking faults, nor regions of misfit as pre- 
dicted for displacement spikes. Since the lattice 
strains around a vacancy or a vacancy cluster 
supposedly are too small to cause a noticeable 
contrast in the electron image, no statement can 
be made as to the presence or absence of these 
lattice defects. Further, the high jump frequency 
of vacancies at pile temperature makes it very 
unlikely that foils 1000 A thick contain vacancies 
after irradiation. 

For these reasons, and in order to test the 
hypothesis derived from the above result that 
the interactions of the first glide dislocations 
with vacancies or vacancy clusters result in 
obstacles to slip, the following experiment was 
conducted. Nickel foils approximately 0.1 mm 
thick were subjected to a neutron flux of 10°* nvt, 
then deformed about 2%, and subsequently pol- 
ished to a thickness of 1000 A and examined. 
Figure 1 is typical for the results obtained. The 
micrograph exhibits tangles of intricately kinked 
dislocations and dislocation loops with diameters 
between 50 A and a few hundred angstroms. This 
may be compared with a corresponding micro- 
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FIG. 1. Electron micrograph of irradiated nickel 
specimen (10'* nut) deformed before thinning to 1000 A 
thickness. Lines represent dislocations, dark spots 
small dislocation loops. 40000:1. 





FIG. 2. Transmission micrograph of deformed non- 
irradiated speciman. 40000:1. 


graph of an unirradiated deformed specimen, 
Fig. 2, which shows a similar dislocation pattern 
but on a coarser scale, and practically no small 
dislocation loops as represented by the dark 
spots in Fig. 1. 

The following interpretation of these findings 
is suggested. The vacancies produced through 
irradiation in 0.1-mm thick specimens form 
larger or smaller vacancy aggregates. On plas- 
tic deformation the first glide dislocations 
moving through the lattice cause the transforma- 
tion of vacancies into condensed layers by a 
mechanism described elsewhere.** These loops 
are nucleated either in the neighborhood of dis- 










locations, or in direct contact with them, in the 
latter case causing irregular kinks in the dis- 
location lines as visualized by Kimura, Maddin, 
and Kuhlmann-Wilsdorf for quenched-in vacan- 
cies in copper.“* In this way numerous con- 
densed loops are formed, distributed throughout: 
the deformed specimen, and many kinked and 
intertwined dislocations. The greatly increased 
yield strength of irradiated nickel specimens at 
room temperature then is not directly caused by 
vacancies and vacancy clusters as present in 

the metal after irradiation but is due to the 
anchoring of the dislocations through loops and 
kinks nucleated at them. After their passage 
along a certain minimum number of slip planes, 
an appreciable portion of the excess vacancies 
will have been eliminated from the regions 
between these slip planes. Thus, paths will have 
been cleared for later dislocations which may 
originate from suitable loops or kinks, and con- 
sequently a low rate of work hardening may be 
expected for irradiated metals, as indeed has 
been observed. 
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SPIN-LEVEL INVERSION AND SPIN- TEMPERATURE MIXING IN RUBY* 
R. H. Hoskins 


—_ 


Research Laboratories, Hughes Aircraft Company, Culver City, California 
(Received June 9, 1959) 


We have performed experiments on the succes- 
sive inversion of spin levels of Cr** in ruby by 
adiabatic fast passage followed by a pulsed mag- 
netic field. The preliminary results reported 
here were all obtained at 1.4°K by using a cavity 
completely filled by solid’: 75% clear sapphire, 
25% ruby.” The pulsed field, in the form of a 
single sinusoid of period 1.6 x10-* sec, was pro- 
duced by discharging a condenser bank into a 
Helmholz pair, 3 cm in diameter, located inside 
the helium Dewar. Peak fields of 18000 oersteds 
could be achieved. The experiments to be dis- 
cussed may be considered to have taken place in 
a time much shorter than 7, ~0.1 sec, the ther- 
mal relaxation time. 

The energy levels appropriate for the crystal 
orientation used in one set of experiments are 
shown in Fig. 1." The experiment is started with 
a static field H, of about 6000 oersteds and the 
magnetic field is pulsed as shown at the bottom 
of Fig. 1. An inverting pulse at 9 kMc/sec serves 
to invert successively the 3-2 and 2-1 transitions* 
as the field is swept down. As the pulsed field 
reverses direction, stimulated emission of radia- 
tion from the 1-2 transition can occur as the field 

















OSCILLATION PULSES 





t=1200 psec 


FIG. 1. Energy levels in ruby near 60-degree 
orientation. The magnetic-field pulse is shown at the 
bottom of the figure. The levels have been numbered 
to agree with the matrix element calculations in refer- 
ence 4; level 1 has the highest energy. 
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passes through the values where the resonant 
frequencies of the cavity are equal to the Larmor 
frequency of the 1-2 transition. Depending upon 
the coupling to the particular cavity mode, either 
oscillation or amplification of a probing signal at 
the appropriate frequency is observed. Amplifica- 
tion has been obtained at frequencies near 14 and 
24 kMc/sec, and oscillation pulses have been ob- 
served near 14, 18, 24, and28 kMc/sec. Peak 
power emitted at 14 kMc/sec was several milli- 
watts. 

The expected rapid attainment of a Boltzmann 
distribution of spin populations when the energy 
levels become equidistant® in a multilevel spin 
system such as Cr*® in ruby has been demon- 
strated directly. Consider the energy-level dia- 
gram in Fig. 2.5 The 9-kMc/sec absorption spec- 
trum obtained with a single sinusoidal field sweep 
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FIG. 2. Energy levels in ruby near 70-degree 
orientation. The transitions shown are 9 kMc/sec. 
The levels have been numbered to agree with the 
matrix element calculations in reference 4; level 1 has 


the highest energy. 
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FIG. 3. 9-kMc/sec absorption spectrum in ruby 
near 70-degree orientation. 


having an amplitude of approximately 2000 oer- 
steds starting at an H of approximately 2400 oer- 
steds, is shown in Fig. 3(a). The field initially 
sweeps upward so that the transitions are tra- 
versed in the order 2-3, 3-4, 3-4, 2-3, 1-2, 
2-3’, 2’-3’, 1-2, the total time being 1.6 milli- 
seconds. In Fig. 3(b) a 9-Mc/sec saturating 

pulse is applied at the time the field first passes 
through the 2-3 transition. The 3-4 absorption 

is enhanced because of the fewer number of spins 
inlevel 3 (points A). The missing absorption at 
point B is the return through the 2-3 transition 
and it shows that this transition is indeed satu- 
rated. The 1-2 transition (point C) is seen to be 
enhanced because of the addition of spins to 

level 2. Nothing untoward has happened thus far. 
However, at points D, the 2’-3’ transition, we 

see that the populations of levels 2 and 3 are no 
longer equal and have, in fact, been restored in 
atime << T, to a distribution comparable to that 
before saturation; furthermore, the 1-2 absorp- 
tion (point E) is no longer enhanced. The reason 
is that the field passed through the point indicated 
by H* in Fig. 2 where the levels 1, 2, and 3 are 
equidistant. In this region, spin-spin interaction 
can set up a Boltzmann distribution between 
levels 1, 2, and 3 in a time of the order of 7, 









which we can estimate from the spin density to 
be less than one microsecond. 

It was found to be impossible to carry out a 
three-step “staircase” inversion scheme in ruby 
near the 70-degree orientation because of this 
mixing. Traversal of enough of the 1-2 transition 
to allow inversion by fast passage brought the 
field into close enough proximity to H* to erase 
the inversion. Thus, in extending the “staircase” 
scheme to more than two steps, difficulty may be 
met not only in cross-relaxation between nearly 
equal pump transitions, ® but also in avoiding field 
regions of equally spaced energy levels. This 
difficulty may be side-stepped by using a pump 
frequency that is high compared with the crystal 
field splitting; in this case, however, a price 
must be paid in the form of smaller transition 
probability for the “forbidden” transition. 

The author would like to thank Dr. G. Birnbaum 
and Dr. T. H. Maiman for helpful discussions and 
C. R. Duncan for technical assistance. 
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FISSION FRAGMENT DAMAGE IN ZIRCONIA 


M. C. Wittels and F. A. Sherrill 
Solid State Division, Oak Ridge National Laboratory,” Oak Ridge, Tennessee 
(Received May 6, 1959) 


Natural monoclinic ZrO, was reported’ to 
transform into the high-temperature cubic modi- 
fication upon exposure to fast neutrons. In that 
same study it was suggested that thermal spikes 
contributed to the transition by means of a rapid 
high-temperature quench. Further efforts to 
study these effects in synthetic single crystals 
and powders of ZrO, have proved unsuccessful 
because the transition could not be effected even 
with irradiation doses higher by a factor of four 
over the original irradiation dosage. This is 
demonstrated in the single-crystal x-ray diffrac - 
tion pattern in Fig. 1. 

Unirradiated portions of the original natural 
material were subjected to neutron activation 
analysis in order to determine the concentration 
of natural uranium and thorium which these 
oxides often contain as impurities. As expected, 
it was found that the material which undergoes a 
phase transition upon fast neutron irradiation 
contains uranium and thorium in concentrations 
of 3550 ppm and 39.2 ppm, respectively. As will 
be shown later, the number of fission events 
which took place in this material upon exposure 
to 9.4x10** neutrons/cm? was of sufficient mag- 
nitude to completely transform the monoclinic 
structure into the high-temperature cubic phase. 
It is apparent, also, that the monoclinic ZrO, is 
stable under fast-neutron bombardment, contrary 
to earlier beliefs,’ and that the observed irradia- 
tion- induced phase transition is produced by the 
action of fission fragments. As a result of these 
findings further experimental evidence is avail- 
able to support proposals of a “fission spike” 





(a) (d) 
FIG. 1. Zero level precession photographs — (HK0) 
projection (a) Before irradiation; (b) After irradiation 
with 3.6 x10" nvt,. 
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mechanism.” ~° 

Samples of ZrO, with known concentrations of 
fissionable impurities were subjected ta thermal] 
neutron dosages at irradiation levels capable of 
producing partial transformation from the mono- 
clinic to the cubic structure. Since these fission- 
able atoms were contained within the crystal lat- 
tice of ZrO, in amounts less than 0.5% and since 
the range of a 100-Mev fission fragment is pro- 
bably less than 8x10~* cm in this material, it is 
reasonable to assume that (1) nearly all of the 
fission fragment energy is dissipated within the 
crystal lattice, and (2) the local areas surround- 
ing each fission event are sufficiently separated 
to be termed independent regions. Then if there 
are no annealing-type rearrangements of atoms 
within the damaged regions a reasonable esti- 
mate of the number of atoms affected by a fission 
fragment can be made by determining the fraction 
of transformed material in partially transformed 
crystals. Figure 2 shows some x-ray diffraction 
patterns from which these estimates are made. 

The method used in this analysis was to meas- 
ure the integrated line intensities of the two 
phases directly from the Debye-Scherrer photo- 
graphs. From the relative change in intensities 
of these lines with increasing irradiation dosage 
the estimates of percentage transformation were 
obtained. It is felt that the over-all accuracy of 
these estimates is within + 25%, the greatest in- 
accuracy being caused by radioactivity - induced 
film fogging. Some overlapping of damaged re- 





(a) (d) (c) (a) (e) 


FIG. 2. Debye-Scherrer photographs of ZrO, irra- 
diated with fission fragments. (a) Before irradiation; 
(b) After irradiation with 7.66 x10 fissions per cm’, 
5 % transformed; (c) After irradiation with 6.67 x108 
fissions per cm’, 50% transformed; (d) After irradia- 
tion with 1.07 x10" fissions per cm’, 80% transformed, 
() After irradiation with 3.3910" {ssions per cm’, 
completely transformed. 
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gions probably occurs with dosages in excess of 
5x10** fissions per cm*, so greater accuracy is 
obtained when less than 50% of the material has 
been transformed. By this analysis it is seen 
that ~ 10° atoms are affected by each fission 
fragment in inducing the phase transition. In the 
original observation of the irradiation- induced 
phase transition,’ the particle size of the trans- 
formed crystallites was ~100A radius as deter- 
mined from x-ray diffraction line breadths. 
Particles of this size would contain about 3.5x10° 
atoms, lower by a factor of three from the 10° 
atoms estimated in the present work. The dis- 
crepancy is not unreasonable, however, since 
the x-ray line breadth determination did not con- 
sider the strain contribution to the line broaden- 
ing. If the true particle size was as much as 
140A radius, no discrepancy would exist and 

both cases would give 10° atoms affected by each 
fission fragment. Other irradiation-induced 
phase transitions in reactor fuel materials®” 
have been reported, and it is interesting to note 
that estimates’ of the number of atoms affected 
ina U-Mo alloy by a single fission event is fairly 
close to the number proposed here for the insula- 
tor ZrO,. Calculations using the method of Seitz® 
show that a pair of fission fragments in a mate- 
rial like ZrO, would produce about 10* displaced 
atoms as a result of elastic collisions. It does 
not seem likely that a single displaced atom 

could affect 200 surrounding atoms to produce a 
phase transition, and this is supported by the 
stability of the material under fast neutron bom- 


bardment, solely. It is more likely that the pro- 
cess under which the transition occurs is a 
“fission spike” mechanism with a rapid, high- 
temperature quench. 

The results of this investigation indicate that 
monoclinic ZrO, is stable under fast neutron 
bombardment but that it transforms into the 
cubic phase under the action of fission fragments. 
In this material it appears that ~2x10° atoms 
are quenched into the high-temperature modifica- 
tion as a result of the energy dissipated by a pair 
of fission fragments. 
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PROTON CAPTURE GAMMA RAYS IN THE GIANT RESONANCE REGION 


H. E. Gove, A. E. Litherland, and R. Batchelor* 
Atomic Energy of Canada Limited, Chalk River, Ontario, Canada 
(Received July 20, 1959) 


Measurements of the 90° yield of gamma rays 
resulting from the capture of protons, with ener- 
gies up to 11 Mev, in B™ and Al” have been made 
using the Chalk River Tandem Accelerator.’ In 
both cases the gamma rays investigated were 
those leading to the ground and first excited states 
(4.43 and 1.78 Mev) of the appropriate residual 
mucleus, C’* and Si**, respectively. In addition, 
the angular distributions of these gamma rays 
with respect to the incident proton beam were 
measured at several proton energies. This Let- 
ter reports principally the results of the reac- 
tion BY*(p, y)C!? (Q =15.949 Mev). 


The gamma rays were detected in a 5-inch 
diameter by 6-inch long unshielded Nal(T1) crys- 
tal.? For the angular distribution measurements 
a second 5-inch diameter by 4-inch long NalI(T1) 
detector was used at a fixed angle as a monitor. 

Both the aluminum and boron targets were thin 
self-supported films through which the beam 
passed to be stopped in a lead plate behind shield- 
ing blocks about 16 feet from the target. The 
boron targets were prepared by mixing finely 
powdered boron of normal isotopic abundance 
with a colloidal graphite suspension in alcohol. 
This mixture was spread on a glass plate and the 
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dried films were floated off in water. The result- the yield of 7, shows a pronounced broad peak 



































ing targets were about 20 kev thick for 10-Mev with a maximum at approximately 22.5-Mev ex- 
protons. Equivalent thicknesses of aluminum citation energy in C’? and a width at half maxi- 
foil were used. mum of about 3 Mev in the center-of-mass sys- 
Figure 1 shows the 90° yield of the ground-state tem. This is consistent with photonuclear giant 
proton capture gamma rays (7,) from the B"(p, resonance observations* and measurements of 
y)C™ reaction as a function both of the incident the B"(p, 7,)C reaction over a limited energy 
proton energy in the laboratory system and of the range by Gemmell et al.® The yield of 7,, on the 
excitation energy in C’*. The points shown rep- other hand, is quite different showing no pro- 
resent several runs with different targets and the nounced maximum in the yield up to the highest 
fluctuations are partly due to the normalization proton energy of 11.4 Mev at which measurements 
required and partly because of target nonuniformity. were made. 
The yield of gamma rays leading to the first ex- Angular distributions of 7, in the B"(p, 7)C” 
cited state of C’’ (these gamma rays will be re- reaction at three energies on the giant resonance 
ferred to hereafter as y,) shows a general rise and one below it are shown in Fig. 2. The three 
between proton energies of 5.5 to 11.4 Mev with measurements on the resonance can be fitted 
some evidence for two or three broad resonances. within the experimental uncertainties (+ 10%) by 
At proton energies corresponding to the peak of an expression of the form W(@) =1 - 0.5P,(cosé) 
the giant resonance for y, the yield of y, is appre- while that at the lower energy is spherically sym- 
ciably lower than that for y,. Consequently a metric. Distributions of y, at the same four 
detailed yield curve of y, awaits the determina- energies are all spherically symmetric within 
tion of the spectrum line shape of high-energy the experimental uncertainties (425%). The ob- 
gamma rays in the 5-inch by 6-inch Nal(Tl) crys- _ served ground-state distribution is that expected 
tal. A careful search, in approximately 20-kev for the case of d-wave spinless particles cap- 
steps, has been made between proton energies tured by a spinless nucleus and then dropping to 
of 5.5 to 7.8 Mev for one of the breaks observed* a p shell with the emission of elect: ‘c dipole 
in the integral yield of neutrons, from brems- radiation. It can also, of course, be accounted 
strahlung irradiation of C’, at 21.72-Mev exci- for by one of the terms in conventional angular 
tation (this corresponds in the present experi- correlation theory. The former interpretation 
ment to a laboratory proton energy of 6.3 Mev). corresponds to direct capture and has been pre- 
This break can be interpreted® as a resonance at dicted for giant-resonance transitions by Courant.’ 
this energy with a total width of 170+ 70 kev. The situation for the Al?"(p, 7)Si”® reaction (Q 
Neither the yield of 7, nor y, showed any varia- =11.588 Mev), which is quite different, will be 
tion from a fairly uniform rise greater than +10% 
over the region in question. As seen from Fig. 1 T T T T T T ——T T 
2 
12 3 " 7 
EXCITATION ENERGY IN C -Mev > 
6 17 18 19 20 2) 22 23 24 25 26 27 a 5 7 
a T T T T T T T T 3 ° 
4 4 _@ ~ oe 
16 B'(p,%)C" 2 
14 4 3 6° - 
S 12+ : = 
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™ 10+ 4 g 
> 1 : 
= 8b : E 
< = - 1 iit L l l l 1 1 1 
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FIG. 2. The angular distribution of gamma rays 
7 leading to the ground state of C” from the B'"(p,y)C” 
fe) SS eS ee ee ee ee ee ee | reaction at a number of proton energies; closed cir- 
’ PROTON ENERGY ee we cles 3.55 Mev, open seuten 6.09 i. crosses 7.21 
Mev, and triangles 8.45 Mev. The solid line is the 
FIG. 1. The 90° yield of gamma rays leading to the function W(@) =1-0.5P,(cos@) and the dashed line is the 
ground state of C from the B'(p,y)C™ reaction. mean of the values for 3.55 Mev. 
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summarized briefly here. The 90° yields of both 
yo and 7, in this case show considerable, though 
dissimilar, sharp resonance structure. The 
envelope of this structure for 7, rises steadily 
for proton energies from 4.0 to 10.4 Mev while 
that for 7, rises and stays relatively constant 
from 6.0 to 10.4 Mev. If, indeed, the ground 
state gamma-ray yield falls off at energies above 
the range presently covered to give a giant reso- 
nance envelope, the present measurements indi- 
cate that the width of the envelope is at least 

4.5 Mev. 

Angular distributions of both 7, and >, in the 
Al?"(p, y)Si*® reaction have been measured at pro- 
ton energies of about 6.9, 7.2, and 8.4 Mev. All 
three distributions for 7, are fitted by W(@) =1 
+a,P,(cos@) where a, is slightly greater than zero, 
while those for 7, are spherically symmetric 


within experimental uncertainties. Direct cap- 
ture processes give rise to a coefficient of P,(cosé 
which is less than or equal to zero and hence the 
observed positive coefficient for 7, suggests that 
the process cannot be described in this way. 





“Seconded from the Atomic Weapons Research Estab- 
lishment, Aldermaston, England. 

‘Designed and built by High Voltage Engineering Cor- 
poration, Burlington, Massachusetts. 

?Manufactured and canned by the Harshaw Chemical 
Company. 

3L.. Katz (private communication). 

“See for references M. E. Toms, “Bibliography of 
Photonuclear Reactions,’ Naval Research Laboratory 
Bibliography No. 14, August, 1958 (unpublished). 

‘Gemmell, Morton, and Titterton, Nuclear Phys. 10, 
33 (1959). 

SE. D. Courant, Phys. Rev. 82, 703 (1951). 





NEUTRON ASYMMETRY FROM MU CAPTURE IN MAGNESIUM 


W. F. Baker and C. Rubbiat 
Columbia University, New York, New York 
(Received July 27, 1959) 


It has been pointed out by many authors’ that 
a measurement of the angular distribution of 
neutrons emitted after the capture of polarized 
u~ mesons by complex nuclei may give infor- 
mation on the nature of the coupling constants in 
the process 


i +p—-n+v. (1) 


A spatial asymmetry with respect to the » spin 
direction would also be evidence for parity non- 
conservation in reaction (1). Until now no such 
violation has been experimentally observed. A 
calculation by Primakoff suggests that an asym- 
metry parameter as large as -0.4 could be ex- 
pected.” In addition to the neutrons from (1), 
evaporation neutrons would be emitted. Since 
these would mask an asymmetry, they were 
discriminated against in the present experiment 
by detecting only those neutrons with an energy 
>5.5 Mev. 

A previous attempt to measure this asymmetry 
was hindered by a large contamination of nuclear 
de-excitation y rays.® A large number of such 
y Tays per p capture has been observed.°™ In 
order to reject these 7 rays in this experiment, 
the neutrons were detected with a sandwich 
counter consisting of 15 6-in. by 3-in. sheets of 
1/8-in. thick plastic scintillator. The adjacent 


layers were viewed alternately by two phototubes 
[6 and 7 of Fig. 1 (a)]. Electrons of energy above 
a certain threshold produced by gamma rays 
crossing the counter have in general a range 
sufficient to penetrate more than one layer thus 
causing a (67) coincidence. It is therefore 
possible to reject such events because most of 
the neutrons from »™ capture produce counts 
either in 6 only or in 7 only, the range of the 
proton recoil being short compared to the thick- 
ness of a single layer (see Fig. 2). 





7 BOLYETHYLENE -. 























TIME FROM MU STOPPING IN TARGET (» SEC) 


FIG. 1. (a) Experimental arrangement. (b) Timing 
diagram. 
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FIG. 2. Proton recoil spectrum from 14.8-Mev 
neutrons in the presence of a high gamma-ray 
background. Upper points are without (67) anticoin- 
cidence requirement; lower points are with this 
requirement. Solid line is calculated proton recoil 
spectrum corrected for loss of protons which penetrate 
more than one layer of neutron counter. 




















The y)-ray rejection efficiency of the counter 
was determined by counting 2.76-Mev (Na™*) and 
4.43-Mev (p+C?~p+C™+y) y rays witha 
pulse-height threshold of 2.2 Mev, with and 
without the anticoincidence (67). The y-ray 
rejection efficiency was found to be 0.96 for 


In our counter, an electron of 2.2 Mev gives the 
same pulse height as a 5.5-Mev proton.® 

The experimental arrangement is shown in 
Fig. 1 (a). A neutron event is defined by a count 
in 6 or 7 in anticoincidence with (a) a count in 5, 
(b) a (67) coincidence, and (c) a (16) or (17) 
coincidence to reject neutrons from pion stars in 
the absorber. The neutron event was delayed 
0.3 usec before triggering two identical 1-ysec 
wide gates. A stopped yu pulse, 234, delayed 
1.8 psec was sent to one gate to determine an 
event, and an undelayed pulse was sent to the 
other to determine an accidental count [Fig. 1 (b)]. 








pulse height of recoil proton. 





2.76-Mev and better than 0.98 for 4.43-Mev » rays. 


Thus both events and accidentals were measured 
simultaneously. The outputs of these gates trig- 
gered two 100-channel pulse-height analyzers on 
which was recorded the total energy lost in the 
neutron counter for events and accidentals. 

The neutron asymmetry was measured by pre- 
cessing the ».~’s stopped in a magnesium target 
through 90° by a constant magnetic field which 
was reversed every 30 minutes. Approximately 
27000 events and 15000 accidental counts for 
each direction of the field were recorded. Only 
accidental counts were observed (a) when the 
target was removed, (b) when a lithium target 
was used, and (c) when y*’s were stopped in the 
target. 

By removing the (67) anticoincidence require- 
ment on the sandwich counter, y rays will be 
counted in addition to neutrons. An increase of 
3.8 times was found in the events rate, corrected 
for accidentals. Combining this with the y-ray 
rejection efficiency above yields a negligible 
contribution to the final result. 

To determine the »~ polarization in the target, 
electrons from yp decays were recorded under 
identical conditions except that a (567) coin- 
cidence was required. The results are summa- 
rized in Table 1 where a is the asymmetry 
parameter in the assumed 1+a cosé@ angular 
distribution. These results seem inconsistent 
with the significantly large negative value that 
has been theoretically predicted. 

We wish to acknowledge the assistance of 
Dr. T. Coffin, Dr. H. Gursky, and Dr. D. Tycko. 
We also want to thank Dr. Lee Lidofsky for the 
use of the Columbia Van de Graaff for calibra- 
ting the neutron counter, and Dr. D. Tycko for 
assistance in analyzing the data. 








Table I. Measured asymmetries in magnesium for neutrons and electrons. 
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We have observed angular correlation effects 
between pions emitted from antiproton annihilation 
events. This experiment was carried out with a 
separated antiproton beam? of momentum a) 
=1,05 Bev/c. A total of 2500 annihilation events 
were observed in 20000 pictures taken with the 
Lawrence Radiation Laboratory 30-in. propane 
bubble chamber. 

Pion pairs formed by the charged pions emitted 
inan antiproton-annihilation event can be con- 
sidered in two groups: viz., like pairs (in the 
isotopic-spin state J=2) and unlike pairs (in the 
isotopic-spin states /=0,1, or 2). We searched 
for correlation effects in these separate groups. 
Our results show that the distribution of the angles 
between pions of like charges is strikingly differ - 
ent from the distribution of the angles between 
pions of unlike charges. The angles between pion 
pairs were computed in the center of mass of the 
antinucleon-nucleon system.” The results shown 
in Fig. 1 were obtained from the analysis of the 
“hydrogenlike” events in which four and six 
charged pions, respectively, are emitted. We 
define as “hydrogenlike” those events giving rise 
oan equal number of positive and negative pions. 
Events showing visible evaporation prongs are ex- 
cluded from this sample. The curves shown in 
Fig. 1 were calculated on the basis of the statisti- 
talmodel, expressed in the Lorentz-invariant 
phase-space* (LIPS) form, for pion production 
from a nucleon-antinucleon annihilation. This 
model imposes energy and momentum conserva- 
— tion, but no other constraints. The distribution 
ofthe pion-pair angles @,, for an annihilation into 
"pions of mass up is* 


$n(COS®,,) = { if P:P2Fy_,(W!)dwydw,, 








| 





vith integration limits from w, 2H, w, >» to max 





PION-PION CORRELATIONS IN ANTIPROTON ANNIHILATION EVENTS* 


Gerson Goldhaber, William B, Fowler, Sulamith Goldhaber, T. F. Hoang, 
Theodore E. Kalogeropoulos, and Wilson M. Powell 
Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received July 17, 1959) 
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FIG. 1. Distribution of angles between pion pairs 
as a function of cos@;,. The curves correspond to 
calculations on the Lorentz-invariant phase-space 
(LIPS) model. The deviations of the experimental dis- 
tribution from the LIPS model are discussed in the 


text. 
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Table I. The ratio y for like and unlike pion pairs and for the Lorentz-invariant phase-space (LIPS) model .2 














Like pions Unlike pions All pions combined Statistical mode] 


No. of pairs Y No. of pairs Y No. of pairs Y Y 





702 1,23+0.11 1404 2.06 0.12 2106 1.72 +0.08 1.74 
214 1.06 +0.15 318 1.91 +0.23 532 1.50+0.13 1.60 








“The ratio y is the number of pion-pair angles greater than 90° compared to those smaller than 90°. The errors 
quoted are the standard deviations based on the number of pairs observed. 





values given by W’” =(n -2)?u”. Here we define 


Ww’? =(W - w, - w.)? - (P, +)’; 


F,,-,(W’”) is the Lorentz invariant phase space (O)LIKE 49 
for (n -2) pions, W is the total energy in the 
center of mass of the antinucleon-nucleon system, 
and we have 





COSO,2 =P,*P2/|P,! | Pe! - 

To compare with the experimental distributions 
for events with n, charged pions, averages over anuenet 4ee 
¢, Values with n >n, are required. This takes 
into account the presence of additional neutral 
pions in the annihilation. We have used the fre- 
quency distribution of the pion multiplicity in 
annihilation events, as reported elsewhere,° for 
computing these averages. 

In Table I we have expressed the distribution 
of pair angles in terms of the ratio, y, of the 
number of pion-pair angles greater than 90° to 
the number smaller than 90°. As can be seen 
from Fig. 1 (c) and (f), the pion-pair distribution 
of like and unlike pions combined agrees very well 
with the LIPS model.® The distribution of angles 
between pions [Fig. 1 (a) and (d)] deviates dis- 
tinctly from the LIPS model. The 74;,. values 
for 4n* and 67+ differ from yy, pg by 5 and 3.4 
standard deviations, respectively, in the direc- 
tion of greater isotropy. The distribution of pion- 
pair angles for unlike pions appears to be slightly 
more asymmetric than the LIPS model predicts. 
In this case, the values of yynjike are 2 and 1.5 
standard deviations, respectively, removed from 
the value given by the LIPS model. 

We have also computed the invariant quantity 
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FIG. 2. The distribution of @, the square of the 
Q12” = (w, + wp)? - (P, +p,)”, total energy in the center of mass of the pion-pion sys 


tem in units of p’. 
for each pion pair. Here Q is the total energy in Te 


the center of mass of the pion-pion system. These 

distributions are given in Fig. 2. Within statisti- been observed. Curves shown in Fig. 2 were als0 

cal limits no significant difference between the computed on the basis of the LIPS statistical model. 
@ distribution of like and unlike pion pairs has The experimental Q’ distributions show no marked 
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deviation from the calculated curves. 

In view of the remarkable agreement of the com- 
bined data with the LIPS model, we consider this 
model to be a good description of the over-all 
— physical situation. We ascribe the deviations 

observed in the like and unlike pion pairs to the 
presence of additional pion-pion correlations. 
— It should be noted that in order to retain the 
over-all agreement with the LIPS model, a cor- 
relation effect present in either of the two distri- 
butions must reflect on the other distribution. 
The effect we have observed clearly indicates the 
need for a refinement of the LIPS model. 

Work is in progress to investigate modifications 
of the LIPS model by means of correlation func- 
tions which may account for the observed pion- 
pion correlations.’ These correlation functions 
involve the radius R of the interaction volume as 
aparameter, and may enable one to determine 
its value. The role of such correlations will be 
to effectively enhance specific states. In this 
connection the recently suggested resonance of 
the 2x system with J=1, J=1 may be of particu- 
lar interest.® 
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In the past year evidence has been accumula- 
ting that high-energy nucleons arrive at the 
earth following solar disturbances.'~* The char- 
-_ ateristics of the arrival of these beams from 
| the sun are the enhanced cosmic noise absorp- 

tion on the polar cap, the arrival of particles at 

the top of the atmosphere, followed by a Forbush 
tecrease in cosmic rays at sea level. 

We report here the most spectacular such 

vent which has occurred to date, in which ion 
tambers, Geiger and scintillation counters, 

ad photographic emulsions were exposed above 


— 





™ Minneapolis at a pressure altitude of 10 g/cm? 
throughout the period of enhanced intensity. We 
iis wreviously reported an event of this type in 


thich the intensity of the integral flux of protons 
ithe top of the atmosphere in a vertical direc- 
lon was increased by about a factor of two, on 
e als0 @ March 26, 1958.? In the present event, the in- 

11 mode '¢gral flux of particles at the top of the atmos- 
marked § there increased by approximately a factor of 










PROTONS FROM THE SUN ON MAY 12, 1959 


E. P. Ney, J. R. Winckler, and P. S. Freier 
University of Minnesota, * Minneapolis, Minnesota 
(Received July 17, 1959) 


1000 above that of cosmic rays, and the composi- 
tion of the incoming beam as observed at 10 g/cm? 
atmospheric depth was essentially pure hydrogen. 
The flux of alpha particles and heavy nuclei was 
not increased and corresponds to the normal 
cosmic-ray flux at solar maximum. However, if 
the solar- injected heavy particles had the same 
rigidity spectrum as the protons we observe, 

the atmospheric cutoff would not allow their de- 
tection at 10 g/cm* atmospheric depth. The 
abundance of relativistic protons and high-energy 
electrons is less than 10% of the flux of solar 
protons. 

The energy spectrum of the incoming particles 
has been measured in photographic emulsions 
and is well represented by a differential kinetic 
energy spectrum of the form n(E)dE =KE~*:* dE 
in the measured energy range 110 Mev<E <220 
Mev. The corresponding integral rigidity spec- 
trum is N(>R) = 7500/R** with R expressed in 
Bv and M(>R) given in protons/m’ sec steradian. 
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This spectrum is quite similar to that inferred 
by Anderson for the August 22, 1958, event® and 
is also like that of the large cosmic-ray flare of 
February 23, 1958.‘ It is, however, much 
steeper than the normal cosmic-ray spectrum 
for which the integral rigidity spectrum is N(>R) 
=K/R?, 

The sequence of events was as follows. At 
2000 UT on the 10th of May, a Class 3+ flare 
occurred at a solar latitude of 19° and a solar 
longitude of 50° in the first quadrant. This 3+ 
flare was the largest flare reported in this solar 
cycle. It lasted for over three hours and gave 
rise to a very strong solar noise storm. About 
a day later, as is usual for these events, there 
was a sudden commencement of a magnetic 
storm which lasted from about 2300 UT on the 
11th of May to 2300 UT on the 12th of May. At 
0030 UT on the 12th of May the Forbush decrease 
on the sea-level neutron monitor at Minneapolis 
began, and reached a maximum of 15% at 0400 
UT on the 12th. The beginning of a Type III 
cosmic noise absorption was observed by H. 
Leinbach at College, Alaska, to occur at about 
0100 UT on the 11th of May. This cosmic noise 
absorption reached a value off scale for Lein- 
bach and exceeding 17 db of absorption through- 
out the latter part of the day on May 11 and up 
until at least 1700 UT on the 12th of May. The 
occurrence of this cosmic noise absorption, 
which is taken now to indicate the arrival of 
particles from the sun, was reported to us by 
Leinbach. We immediately began making balloon 
flights, and in the interval between 0200 UT on 
the 12th of May and 1500 UT on the 12th of May, 
five balloon flights were launched. Of these 
flights, four were recovered, and allow the study 
of the time variation of the particle intensity in 
the emulsions as well as the detailed structure 
of the time variation of the event in the ion cham- 
bers and single counters. The ion chamber 
shows that the ionization at the top of the atmos- 
phere increased to approximately 160 times that 
due to normal cosmic rays. Examination of the 
emulsions showed that the angular distribution 
of the particles at 10 g/cm’ residual pressure is 
very strongly peaked in the vertical direction 
with more than half of the tracks coming from 
within 20° of the zenith. Because of this strongly 
peaked vertical distribution, the vertical flux 
increased by much more than did the total ioniza- 
tion, and the increase in this vertical flux is ap- 
proximately 1000 times cosmic-ray flux, giving 
an integral flux of particles at the peak of this 
event of one million particles per square meter 
184 








per second per steradian, with rigidities greater 
than 0.49 Bv. The angular distribution observed 
in the plates, together with the energy spectrum, 
allow one to determine the angular distribution 
above the atmosphere, and this has been found 
to be isotropic. We therefore feel that it is un- 
likely that the particles observed in this event 
were released from the Van Allen belt, and be- 
lieve that they arrived at the earth as a conse- 
quence of acceleration associated with the 3+ 
flare on the 10th of May. 

It is important to point out that the particles in 
the event, as in the March 26 event, arrived at 
Minnesota at energies below the Stormer cutoff 
energy established for this latitude. It is there- 
fore necessary, in order for us to see them in 
Minneapolis, for the earth’s magnetic field to be 
distorted by the incoming solar beam. The fact 
that the Forbush decrease had reached its maxi- 
mum value at the time that the particles arrived 
at Minneapolis is taken to indicate the arrival of 
the solar.magnetic cloud, which is presumably 
the cause of the lowered magnetic cutoffs which 
allowed these normally forbidden particles to be 
observed in Minneapolis. The arrival of the 
corpuscular beam at Minneapolis was accompa- 
nied by a very strong aurora, and the instrumen- 
tation on the balloon flights was adequate to show 
the presence of the auroral x-rays and to sepa- 
rate these from the nucleonic component. 

The density of ending particles in the emul- 
sions of the four recovered flights has been used 
to study the time variation of the slow proton 
flux. The flux decreases rapidly with time after 
the abrupt beginning of the event, and the inten- 
sity drops a factor of ten every four hours and 
seems to be dying off approximately exponen- 
tially with time. 

Figure 1 shows the counting rate of the various 
counting instruments carried on this flight, as a 
function of Universal Time on the 12th of May, 
together with normalization curves for identical 
instruments flown on the 7th of May at a corre- 
sponding altitude. Figure 2 shows a plot of the 
particle fluxes in this event, together with those 
in the March 26 event, and the rigidity spectrum 
of cosmic rays. The data in Fig. 2 are obtained 
by directly measuring the fluxes of particles at 
various depths in the photographic emulsions. 

This is the most spectacular nucleonic event 
we have observed in two years’ monitoring of 
cosmic rays during the IGY, in which balloons 
were kept at the top of the atmosphere approxi- 
mately 5% of the time. The integral particle 
fluxes observed in this event are higher than 
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counter, and integra- 
ting ionization chamber 
on balloon flight IGC-8. 
For comparison, note 
same instruments on 
flight IGC-7, five days 
earlier. Solar parti - 
cles are first detectable 
at 0430 UT at a balloon 
depth of 100 g/em*. The 
balloon reached ceiling 
at 0530 UT. Note large 
selective response of 
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those reported by Freden and White® for the 

inner radiation belt, but the energy spectrum of 
this event is much steeper than the energy spec- 
trum of the trapped protons in the Van Allen 

belt. 

We wish to express our very great appreciation 
to Harold Leinbach for his prompt warning that 
this event was occurring. It would have been 
impossible to carry out the balloon flights and to 
prepare the equipment without the excellent col- 
laboration of Mr. William F. Huch and Mr. Ray- 
mond W. Maas. 
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FIG. 2. Integral rigidity spectrum of protons ob- 
served during the solar flares of May 12, 1959 and 
March 26, 1958. The cosmic-ray proton spectrum 
‘uring solar maximum and solar minimum is that 
neasured by McDonald [F. B. McDonald and W. R. 
a (to be published); F. B. McDonald (to be pub- 
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SEARCH FOR LEPTONIC DECAYS OF THE SIGMA* 


Jack Leitner,? Paul Nordin, Jr., Arthur H. Rosenfeld, Frank T. Solmitz, and Robert D. Tripp 
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(Received July 17, 1959) 


As part of a continuing study of 8000 K absorp- 
tions in hydrogen’ we have examined ~ 2000 © 
decays for evidence of leptonic modes. 

In order to obtain events which were both unam- 
biguous and accurately measurable, it was re- 
quired that the length of the = track be greater 
than 0.5 mm, and the length of the decay-product 
track be greater than 5 cm and have a dip less 
than 60°. Of the total sample, 750 =~ and 250 
=+~—2++n (effectively 500 £*) satisfied the above 
criteria. 

The search for leptonic decay modes was car- 
ried out by studying the (lab) momentum distri- 
bution of the =-decay products (shown in Fig. 1).? 
The vast majority of the events lies within the 
region ~160 to 224 Mev/c—the kinematic limits 
of the normal £-decay modes, 1*++n. It is there- 
fore clear that careful consideration must be 
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FIG. 1. Histogram of the = decay secondaries. 
The solid curve is the phase-space spectrum of 
= —-e~+n+v. The broken curve is the theoretical 
spectrum of £*—-1x*+n, taking account of experimental 
resolution. 
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given to the problem of momentum measurement ' 
before attempting to estimate the very small lep- 
tonic-mode contribution. 

In order to investigate the effects of multiple * 
and (unobservably) small-angle single scattering 
we studied the momentum spectrum from the re- 
action K- +p- 2 +m* at rest, which yields a 
unique pion momentum of 174 Mev/c. One should 
expect the measured distribution to be symmetri- 
cally distributed about 174 Mev/c (when plotted ! 
against curvature), the spread arising mostly 
from Coulomb scattering of the pion. Multiple 
scattering should produce a Gaussian distribution 
of half-width o, which, for tracks of 15 cm mean 
length, 30° average dip, and 6 =0.8, is o=:11 
Mev/c. The systematics of single scattering 
leads one to expect that about 8 % of the total nun- 
ber of events should have momenta p, outside the | 
interval 174420 Mev/c with a distribution (174 
-p_)*p,.* The calculated momentum spread due 
to the combined effects of multiple and single 
scattering is shown in the dashed curve of Fig. 2? oa 
(plotted against momentum and thus asymmetric). an 
The measured momentum distribution g(p,) of 
these selected primary pions is shown as the 
histogram of Fig. 2. On the basis of the excellen| Firs! 
agreement between the observed and theoretical | dy 
momentum spread, we conclude that our momen-  ‘igib 
tum errors are adequately understood. » mom 

It is convenient to use gb,) as the composite i able. 
error function for momentum measurement. With) mom 
appropriate adjustment to account for the differ- | 5.5 ti 
ence in mean momentum, folding g(p,) into the | We 
known c.m.-lab transformation‘ for the reaction | mom 
=+~1*+n, we obtain the dashed curve of Fig. 1. 5 phase 
The agreement between the observed spectrum | Vé ob 
and that expected from the normal pion decay is | decay 
excellent, indicating that the contribution of other ‘(750 
decay modes is indeed small.° the re 

In the UFI (universal Fermi interaction) theory} o Tal 
of Feynman and Gell-Mann,® leptonic-mode rates) If th 
are ~2 to 6%;” leptons should be emitted witha | bined 


of events 


Number 


phase-space momentum spectrum (shown as the rates 
solid curve of Fig. 1). Since the leptonic spec- a 
e ch 


trum extends well into the pionic spectrum, 4 
detection cutoff must be chosen in order to ob- | Table 
tain upper limits on the rate of leptonic decay. 

The choice of cutoff is somewhat arbitrary; we 
choose it at p, =100 Mev/c for several reasons. 
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9 $ FIG. 2. Histogram of the 7* momentum distribution 
6. | from K +p-~=+n*. The dashed curve represents 


‘_ the calculated spread due to Coulomb scattering. 
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slien First, the background of single-scattered pions 

ical | and pions from radiative decay® should be neg- 

men- ligible below 100 Mev/c. Second, electrons with 

, momenta below the cutoff should be distinguish - 

ite | able (with high efficiency) from pions of the same 
wo momentum, which have an ionization density of 

fer- 3.5 times minimum. 

he | Wehave seen no £-decay secondaries with 

tion momenta below 115 Mev/c. Since the fraction of 

g. 1.5 phase space below the 100-Mev/c cutoff is 33 %, 

ym | We obtain an upper limit for the 2" -electronic 

y is | decay rate, for example, of f([" -e” +n+v) 

other <(750x0.33)-2=0.4%. Branching fractions for 

the remaining modes® are given in the first column 

eory) of Table I. 

rates) Ifthe results of the present experiment are com- 

ha | dined with the world total of published data, the 

the | tates shown in the second column of Table I are 

ec- | dbtained.® In spite of the ambiguity inherent in 

a | te choice of detection cutoff it is clear from 

yb- | Table I that the observed rates are less than 

xy. } ose expected on the basis of the UFI theory by 
we } out an order of magnitude. It has recently been 











ons. | own that a similar discrepancy exists in A° 












Table I. Upper limits to the branching fractions for 
leptonic decays. 








Maximum Rates (%) 


Leptonic mode This experiment Total world data 





Le +ntv 0.4 0.2 
rte +n+v 0.6 0.4 
L>—p-+n+ v 0.4 0.2 
Lt —ptt+nt v 0.5 0.3 
p-—are{e hs y 0.1 

y 
pt —atefehe v 0.2 





leptonic decay.’° While these results in no way 
detract from the enormous success of the Uni- 
versal Fermi theory in nonstrange and K-particle 
weak interactions, it is clear that a simple ex- 
tension of the Fermi interaction with the inclu- 
sion of 27 and An terms cannot be considered an 
adequate model for the description of hyperon 
decays. 





"This work supported in part by the U.S. Atomic 
Energy Commission and the Office of Scientific 
Research. 
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‘Alvarez, Bradner, Gow, Rosenfeld, Solmitz, and 
Tripp, Nuovo cimento 5, 1026 (1956). 

?The data shown in Fig. 1 have been corrected for 
a small in-flight contamination of ~ 3% 

3alvarez, Leitner, Rosenfeld, and Solmitz, Berkeley 
Engineering Report 4310-03 (unpublished) . 

‘Leitner, Nordin, Rosenfeld, Solmitz, and Tripp, 
University of California Radiation Laboratory Report 
UCRL-8737, May 1959 (unpublished). 

5s. Barshay and R. E. Behrends, Brookhaven 
Report BNL-3956, Dec. 1958 (unpublished). 

SR. P. Feynman and M. Gell-Mann, Phys. Rev. 109, 
193 (1958). 

"In the original presentation of reference 6, the 1 
leptonic decay is absolutely forbidden. However, if 
the As=1 selection rule is relinquished, the 4 
leptonic decay is allowed and should occur with a 
relative frequency of ~ 3%. 

®Note that, to within a factor of ~ 2, these upper 
limits depend upon the choice of detection cutoff. For 
example, if the cutoff is taken at 135 Mev/c, about 2 
standard deviations from the minimum momentum 
available to the pion from = —~21*+n decay, the upper 
limit on 2 —e +n+v decay becomes (1 +4)/750 <0 .66 
=0.2%+1.0%. 
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®*The efficiency for the detection of the modes 
x*—A°+e*+ v is 100% because of the low maximum 
electron momentum, leading to the lower rates 
shown in Table I. However, mere phase-space 
considerations indicate an expected rate of ~ 10. 


Crawford, Cresti, Good, Kalbfleisch, Stevenson, 
and Ticho, Phys. Rev. Letters 1, 377 (1958); Orear, 
Nordin, Rosenfeld, Solmitz, and Tripp, Phys. Rey. 
Letters 1, 380 (1958); F. Eisler et al., Nevis Report 
No. 67, March 1958 (unpublished). 





CAUSE OF THE MINIMUM IN THE EARTH’S RADIATION BELT* 


S. F. Singer 
Department of Physics, University of Maryland, College Park, Maryland 
(Received June 29, 1959) 


In earlier publications’ we have put forth the 
hypothesis that the trapped radiation encountered 
by satellites” at low altitudes near the equator 
consists of fast protons; they arise from the 
decay of upward-moving cosmic-ray albedo neu- 
trons which originate in the earth’s atmosphere. 
We also calculated the expected altitude depend- 
ence and energy spectrum of these protons. 

It is our view that solar corpuscular radiation, 
connected with magnetic storms and aurora, must 
exist in the trapped condition at higher latitudes 
and altitudes® and is therefore distinct in origin 
as well as in location from the radiation belt 
found at the equator at low altitudes. 

One of the attractive features of the cosmic- 
ray neutron radiation belt is the fact that its 
injection rate can be accurately calculated from 
cosmic-ray data. Therefore there exists the 
possibility of using the measured intensity to 
determine the density of the outer atmosphere 
of the earth. If only the density of the neutral 
component of the exosphere were to control the 
particle lifetime and therefore the equilibrium in- 
tensity, then the radiation belt intensity should 
increase monotonically with altitude. If, however, 
there are additional factors limiting the lifetime, 
then the intensity would fall off beyond a certain 
altitude and therefore show a maximum, which 
we tentatively placed at 1 to 2 earth radii.’ 

Experimental results from the space probes*»® 
seem to have borne out this hypothesis. Two dis- 
tinct maxima are found: the one at 1.52, (earth 
radii) can presumably be identified with the cos- 
mic-ray radiation belt and the one at 3-4 earth 
radii with the solar radiation belt. 

However, it has been argued that both belts 
may be of solar origin? with the minimum (“slot”) 
at ~2Rp due to some other cause, for example 
an instability of the earth’s magnetic field at 
this location,® or as another possibility the effect 
of the large magnetic anomaly at Capetown.” If 
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either of these views is correct, then we will 
have little chance of deriving atmospheric densi- 
ties from measured radiation belt intensities, 

In order to decide among the various possibili- 
ties, it seemed worth while to calculate the ex- 
pected altitude dependence of the cosmic-ray belt 
by incorporating in more detail the breakdown 
of the adiabatic condition. 

We use the Alfvén discriminant® x =p gradB/B, 
where p is the particle’s radius of curvature, We 
can write x also as (m8c/B)(3M/*)(Mr7*) *«gr, 
Here 8 =V/c and M is the earth’s magnetic dipole 
moment. 

In our calculation we will assume that the life- 
time (and therefore intensity) is primarily de- 
termined by the atmospheric density for particles 
of small x. In other words, at very low altitudes 
the intensity will be primarily controlled by at- 
mospheric density but at higher altitudes, or for 
large momenta, the lifetime is limited, as the 
particle’s pitch angle changes slightly on each 
of its many reflections from the mirror points. 
Just at what point this breakdown of the adiabatic 
condition occurs cannot be easily decided from 
a priori considerations. However, the work of 
Northrop et al. indicates® that the breakdown is 
quite sharp and we therefore take the lifetime to 
be limited by a factor exp[-(x?/b?)]. This func- 
tion drops rapidly when x>b. We can therefore 
derive a maximum value of § at each altitude: 


x7 /”. ) 


Bn 
Bmin 18 of course given by the conditions of the 
experiment and depends on such details as the 
thickness of the counter wall, etc. 
To incorporate this refinement into the earlier 
theoretical discussion” we define the differential 
directional concentration: 


n(B,a,,7,)=H(C,/2n)(Cyy)"8-7*?, 
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where 
C,=9 x, C, =2.6 x10~(Z/A)cd; 


where d(r,, @,) refers to the mean value of at- 
mospheric density experienced by the particle 
during one complete back and forth trajectory 
(the subscript e refers to respective values at 
the equatorial plane); » =2.6; and 7 is an injec- 
tion coefficient™ and of the order of 3. 

We can now calculate as before the integral 
omnidirectional intensity 


9 -14 “1 
1=6.8x10"“fjcosa d “48 Bn) —(3) 


min 
where we take 7)~0.4, cosa, ~0.34, d,(r,) =0.7d, 
but where f now depends on 7: 

f(B ax’ B ain) = bi -4y 2.8 a (4) 

If, e.g., at sea level we take 8,.,, to be 0.7,” 

then f becomes 0.5 at sea level and progressively 
smaller at higher altitudes. If we choose Bmin 

to be 0.2 then f becomes Zero at 1.87 earth radii 
from the center of the earth. 

Using now in Eq. (3) the relative distribution 
with altitude of neutral hydrogen in the earth’s 
exosphere’® (see Table I), we can construct the 
intensity variation with altitude shown in Fig. 1. 
Without any further assumptions the intensity 
shows a maximum at 1.5 earth radii which cor- 
responds closely to the observations. 

The hypothesis that the adiabatic condition breaks 
down at higher altitude can be easily verified. It 
should lead to a progressive diminution of Bax 
inaccordance with Eq. (1) and as seen from Fig. 1. 
This would lead to a progressive “softening” of 
the radiation as one moves out to higher altitudes. 
For example, at 1.5 earth radii there should be 
no protons present with energies greater than 
about 50 Mev and at 1.8 earth radii no protons 
with energy greater than about 25 Mev. (This is 
assuming that 6 is 0.7; obviously higher value of 
bwould lead to higher maximum energies at var - 


lous altitudes.) This point can be checked experi- 
mentally. 


Table I, Altitude variation of Bmax’ J» relative density 
equator. 
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FIG. 1. Omnidirectional integral intensity (relative) 
sf cosmic-ray radiation belt near equator. £,,,, varies 
as indicated; 6, ;,, is taken as 0.2 (20-Mev protons). 
The dashed line Rows the intensity variation expected 
if only hydrogen existed in the exosphere. The pres- 
ence of oxygen, however, causes the steep initial rise 
shown. The transition from O to H, we calculate, 
should take place at about 1000 km. However, it will 
be interesting to obtain an experimental verification of 
this. 


It might be mentioned parenthetically that if this 
theory for the “slot” is correct, then the electrons 
released in high-altitude nuclear explosions (Argus 
tests), which extended out as far as 2 earth radii 
in the equatorial plane, are well within the adia- 
batic limit; therefore their lifetime is controlled 
by atmospheric density rather than by a break- 
down of the adiabatic condition. 
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of exosphere, and relative radiation belt intensity at the 














r/Rp 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 
Bmax 0.7 0.58 0.485 9.415 0.357 0.31 0.273 0.242 0.216 
f 0.5 0.36 0.26 0.19 0.134 0.093 0.061 0.035 0.015 
dy 75 34 23 15 10.5 - 6.1 4.8 
fla, 48 76.5 82.5 89.3 88.5 76 57.5 31 
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®*The efficiency for the detection of the modes 
z*—A°+e*+ v is 100% because of the low maximum 
electron momentum, leading to the lower rates 
shown in Table I. However, mere phase-space 
considerations indicate an expected rate of ~ 10. 
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and Ticho, Phys. Rev. Letters 1, 377 (1958); Orear, 
Nordin, Rosenfeld, Solmitz, and Tripp, Phys. Rey, 
Letters 1, 380 (1958); F. Eisler et al., Nevis Report 
No. 67, March 1958 (unpublished). 
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S. F. Singer 
Department of Physics, University of Maryland, College Park, Maryland 
(Received June 29, 1959) 


In earlier publications’ we have put forth the 
hypothesis that the trapped radiation encountered 
by satellites? at low altitudes near the equator 
consists of fast protons; they arise from the 
decay of upward-moving cosmic-ray albedo neu- 
trons which originate in the earth’s atmosphere. 
We also calculated the expected altitude depend- 
ence and energy spectrum of these protons. 

It is our view that solar corpuscular radiation, 
connected with magnetic storms and aurora, must 
exist in the trapped condition at higher latitudes 
and altitudes® and is therefore distinct in origin 
as well as in location from the radiation belt 
found at the equator at low altitudes. 

One of the attractive features of the cosmic- 
ray neutron radiation belt is the fact that its 
injection rate can be accurately calculated from 
cosmic-ray data. Therefore there exists the 
possibility of using the measured intensity to 
determine the density of the outer atmosphere 
of the earth. If only the density of the neutral 
component of the exosphere were to control the 
particle lifetime and therefore the equilibrium in- 
tensity, then the radiation belt intensity should 
increase monotonically with altitude. If, however, 
there are additional factors limiting the lifetime, 
then the intensity would fall off beyond a certain 
altitude and therefore show a maximum, which 
we tentatively placed at 1 to 2 earth radii.’ 

Experimental results from the space probes*® 
seem to have borne out this hypothesis. Two dis- 
tinct maxima are found: the one at 1.5X, (earth 
radii) can presumably be identified with the cos- 
mic-ray radiation belt and the one at 3-4 earth 
radii with the solar radiation belt. 

However, it has been argued that both belts 
may be of solar origin? with the minimum (“slot”) 
at ~2R, due to some other cause, for example 
an instability of the earth’s magnetic field at 
this location,® or as another possibility the effect 
of the large magnetic anomaly at Capetown.” If 
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either of these views is correct, then we will 
have little chance of deriving atmospheric densi- 
ties from measured radiation belt intensities, 

In order to decide among the various possibili- 
ties, it seemed worth while to calculate the ex- 
pected altitude dependence of the cosmic-ray belt 
by incorporating in more detail the breakdown 
of the adiabatic condition. 

We use the Alfvén discriminant® x =p gradB/s, 
where p is the particle’s radius of curvature. We 
can write x also as (m8c/B)(3M/%*)(Mr-*) *agr, 
Here 8 =V/c and M is the earth’s magnetic dipole 
moment. 

In our calculation we will assume that the life- 
time (and therefore intensity) is primarily de- 
termined by the atmospheric density for particles 
of small x. In other words, at very low altitudes 
the intensity will be primarily controlled by at- 
mospheric density but at higher altitudes, or for 
large momenta, the lifetime is limited, as the 
particle’s pitch angle changes slightly on each 
of its many reflections from the mirror points. 
Just at what point this breakdown of the adiabatic 
condition occurs cannot be easily decided from 
a priori considerations. However, the work of 
Northrop et al. indicates® that the breakdown is 
quite sharp and we therefore take the lifetime to 
be limited by a factor exp[-(x?/b?)]. This func- 
tion drops rapidly when x>b. We can therefore 
derive a maximum value of § at each altitude: 


—_ = wai (1) 
Bmin 18 of course given by the conditions of the 
experiment and depends on such details as the 
thickness of the counter wall, etc. 
To incorporate this refinement into the earlier 
theoretical discussion” we define the differential 
directional concentration: 


n(8,a,,7,)=(C,/2n)(Cyy) 13-747, 
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where 
C,=9x10 —. C, =2.6x10 “3(Z /A)cd; 


where d(r,, ao) refers to the mean value of at- 
mospheric density experienced by the particle 
during one complete back and forth trajectory 
(the subscript e refers to respective values at 
the equatorial plane); » =2.6; and 7j is an injec- 
tion coefficient™ and of the order of 3. 

We can now calculate as before the integral 
omnidirectional intensity 


_ -14 — 
1=6.8 x10 Hjcosa da, ee 


), (3) 


where we take 7j~0.4, cosa, ~0.34, d,(r,) =0.7d, 
but where f now depends on 7: 
f(g ) = b}-47-2.8 -B 1.4 (4) 


max’ B min min 


If, e.g., at sea level we take 8,.,, to be 0.7,” 
then f becomes 0.5 at sea level and progressively 
smaller at higher altitudes. If we choose Bmin 

to be 0.2 then f becomes Zero at 1.87 earth radii 
from the center of the earth. 

Using now in Eq. (3) the relative distribution 
with altitude of neutral hydrogen in the earth’s 
exosphere*® (see Table I), we can construct the 
intensity variation with altitude shown in Fig. 1. 
Without any further assumptions the intensity 
shows 2 maximum at 1.5 earth radii which cor- 
responds closely to the observations. 

The hypothesis that the adiabatic condition breaks 
down at higher altitude can be easily verified. It 
should lead to a progressive diminution of Binax 
inaccordance with Eq. (1) and as seen from Fig. 1. 
This would lead to a progressive “softening” of 
the radiation as one moves out to higher altitudes. 
For example, at 1.5 earth radii there should be 
no protons present with energies greater than 
about 50 Mev and at 1.8 earth radii no protons 
with energy greater than about 25 Mev. (This is 
assuming that b is 0.7; obviously higher value of 
bwould lead to higher maximum energies at var - 
ious altitudes.) This point can be checked experi- 
mentally. 
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FIG. 1. Omnidirectional integral intensity (relative) 
<f cosmic-ray radiation belt near equator. Bmax Varies 
as indicated; 6,;, is taken as 0.2 (20-Mev protons). 
The dashed line dhews the intensity variation expected 
if only hydrogen existed in the exosphere. The pres- 
ence of oxygen, however, causes the steep initial rise 
shown. The transition from O to H, we calculate, 
should take place at about 1000 km. However, it will 
be interesting to obtain an experimental verification of 
this. 





It might be mentioned parenthetically that if this 
theory for the “slot” is correct, then the electrons 
released in high-altitude nuclear explosions (Argus 
tests), which extended out as far as 2 earth radii 
in the equatorial plane, are well within the adia- 
batic limit; therefore their lifetime is controlled 
by atmospheric density rather than by a break- 
down of the adiabatic condition. 





"This research was supported in part by the Geo- 
physics Research Directorate of the Air Force Cam- 
bridge Research Center under contract. 

‘1g. F. Singer, Phys. Rev. Letters 1, 171, 181 (1959). 

*van Allen, McIlwain, and Ludwig, J. Geophys. Re- 
search 64, 271 (1959). 

3g. F. Singer, Bull. Am. Phys. Soc. 1, 229 (1956); 
Trans. Am. Geophys. Union 38, 175 (1957). 

‘J. A. Van Allen and L. A. Frank, Nature 183, 430 
(1959). 

5§. N. Vernov et al., Doklady Akad. Nauk S.S.S.R. 
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"a. J. Dessler, Science 129, 1012 (1959). 

*H. Alfvén, Cosmical Electrodynamics (Oxford Uni- 
versity Press, London, 1950), first edition, p. 33. 

*R. Post, University of California Radiation Labora- 
tory Report UCRL-5044, June 27, 1958 (unpublished), 





p. 6. 
1%» Eq. (10) of reference 1, y should appear as y~'. 
Subsequent expressions must then be reduced by a 


————__ 


factor (2.6)*. 

114 detailed discussion of the injection factor will be 
given elsewhere; it allows one to obtain the latitude 
and angular distribution of the trapped radiation which 
does not concern us in this present note. 

“We will use a nonrelativistic approximation through- 
out; if larger values of g are considered, our treatment 
should be suitably extended. 

13g. J. Opik and S. F. Singer, Phys. Fluids (to be 
published) . 





PARITY OF THE SECOND PION RESONANCE* 


L. F. Landovitz and L. MarshallT 
Brookhaven National Laboratory, Upton, New York 
(Received July 23, 1959) 


We would like to point out in this note an ambig- 
uity in the interpretation of the experiment by 
Stein’ on the polarization of the recoil proton in 
photoproduction in the vicinity of the second reso- 
nance: He found that at 90° in the center-of-mass 
system the polarization was quite large —in fact, 
the polarization seemed to increase to what ap- 
pears to be a maximum at or near the second 
resonance. Following the work of Sakurai,’ it was 
assumed that the polarization is the result of in- 
terference between the first and second resonance, 
and consequently that the second resonance has 
parity opposite to that of the first resonance. 

However, there exists a third resonance** in 
the 7” p system and this should be reflected in 
photoproduction.® The large polarization at 90° 
merely indicates that if the polarization is due to 
the interference of two levels, then these must 
have opposite parity (if they were to have the 
same parity, the polarization would be zero at 
90°). But the interference may be between the 
second and third states with relatively little 
contribution from the “classical” resonance (or 
none at all if the first and second resonances 
have the same spin and parity but differ in iso- 
topic spin). Further, the second and third reso- 
nances are slightly closer to one another than the 
first and second, and the third is broader than 
the “classical” resonance. 

For example, the angular distribution and 
polarization data can be fitted by the following 
level scheme: 


Photon 
x. B. 

Mev T J P 
300 . 32 32 + 
700 12 32 + 

1100 12 32 - 


At the second resonance the angular distribution 
is essentially 1 +(3/2)sin?@ and polarization has a 


190 


maximum at 90° in the center-of-mass system. 
The maximum value of polarization, as pointed 
out by Peierls,® can be 80%. 

The existence of cos*@ and cos*é@ terms (which 
in any case are not large), can be explained by 
the existence of resonances at energies above 
1100 Mev with J=5/2. Alternatively assignment 
of the 1100-Mev resonance to a D,,, level also 
can explain the observed polarization since the 
maximum polarization can be 80% at 90°. 

A caveat to the above alternative possibilities 
is that they cannot explain the sign of the cos’6 
term in 7* production in the region ~500 Mev 
where it has a small positive value.* However, 
this sign persists only in a short energy interval 
and may be spurious. On the other hand, if this 
effect is real, it constitutes evidence for the D,, 
assignment of the second resonance on the basis 
of the assumptions previously indicated, i.e., 
only two states are significant. 

In any case, it seems to us premature to make 
any firm parity assignment to the intermediate 
resonance without further investigation. 





‘Work performed under the auspices of the U. S. 
Atomic Energy Commission. 

Ton leave from University of Chicago, Chicago, 
Illinois. 

‘Pp. C. Stein, Phys. Rev. Letters 2, 473 (1959). 

2J. J. Sakurai, Phys. Rev. Letters 1, 258 (1958). 

3H. C. Burrowes et al., Phys. Rev. Letters 2, 119 
(1959). 

‘Crittenden, Scandrett, Shephard, and Walker, Phys. 
Rev. Letters 2, 121 (1959). 

5F. P. Dixon and R. L. Walker, Phys. Rev. Letters 


1, 142 (1958). 


SF. P. Dixon and R. L. Walker, Phys. Rev. Letters 
1, 458 (1958); R. E. Peierls, thesis, Cornell Univer- 
sity, 1959 (unpublished); F. Turkot, thesis, Cornell 
University, 1959 (unpublished). 
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MODIFICATION OF THE EFFECTIVE-RANGE FORMULA FOR NUCLEON-NUCLEON SCATTERING* 


H. Pierre Noyes and David Y. Wong 
Lawrence Radiation Laboratory, University of California, Berkeley and Livermore, California 
(Received July 27, 1959) 





The effective-range formula is well known to 
be a useful description of low-energy scattering 
phase shifts. However, except at the limit of 
zero kinetic energy, the effective-range formula 
is, in general, not exact. We report here a study 
of the deviation from the effective-range formula 
in light of the analytic structure of partial-wave 
amplitudes as suggested by Mandelstam.’ The 
problem under consideration is the s-wave nuc- 
leon-nucleon scattering in both the singlet and 
the triplet states. 

The effective-range formula, as it stands, 
implies that the only singularities of the s-wave 
amplitude are two poles in the complex momen- 





tum plane: 

1 id , a (q +92) 

1 10 .ind =1/|-+42¢ - iq) =- , 
, \-4-5¢ i) i(q - 9,)(4 - 42) ” 
where 


anft-(-2)), a-ff-(8-2)] 0 


Here g is the center-of-mass momentum, a is 
the scattering length, and 7 is the effective range. 
The effective-range expression (1) can equiva- 
lently be characterized by q, and the residue at 


q: 
Es i ULth 
4, - 2 


where g, is chosen rather than gq, because it has 
acloser relation to the “interaction.” In fact, if 
(i) is sufficiently large, g, will become a bound- 
state pole. This connection between the triplet 
s-wave scattering amplitude and the deuteron 

bound state is well known. We note also that, as 
long as the effective range is positive, q, is al- 
ways on the upper half-plane. 

According to Mandelstam’s representation, 
fartial-wave amplitudes for nucleon-nucleon 
scattering are analytic on the upper-half g plane 
except for branch points on the positive imaginary 
axis corresponding to thresholds for one-, two-, 
ind three-meson exchange, etc. Bound-state 
bles will appear as the “strength” of the branch 
pints becomes sufficiently great. In view of the 
ualyticity implied by Mandelstam’s representa- 
tio, the “interaction pole” (q,) in the effective- 
‘ange formula can be considered as an approxi- 


(3) 





mate replacement for the branch cuts predicted 
by field theory. It is evident that an improvement 
over the effective-range formula can be obtained 
by including the one-meson branch cut exactly, 
and allowing the pole to represent only the aver- 
age contribution of the remaining cuts which are 
farther away from the physical region (the one-, 
two-, three-, ----meson branch points are 
located at q=i/2, i, 3i/2, -++, respectively). 
We shall construct such a function in the follow- 
ing paragraph. 

It is convenient, at this point, to introduce the 
momentum-square variable, v=q*’. The s-wave 
amplitude can be written as 





1 . 
WY) = Taeot6(y) -a™ 


The unitarity condition implies that the inverse 
function h~*(v) has a branch point at y=0 with a 
discontinuity across the cut from 0 to © given by 
- 2ivv.~ The one-meson cut for h(v) can be caicu- 
lated exactly in terms of the renormalized pion- 
nucleon coupling constant f?. The discontinuity 
across the one-meson cut is simply (1if?M/2v) 
for -«<py<-4, M is the nucleon mass in pion 
units. This same cut holds for both the singlet 
and the triplet s-wave amplitudes. The mixing of 
d wave in the triplet state is neglected in this cal- 
culation. An approximate calculation of the s-d 
mixing has already been reported by one of us 
(DYW).? For the construction of a function with 

a given branch cut along the negative real axis 
and a given branch cut for the inverse function 
along the positive real axis, we express h(v) in 
the form of a quotient,° 


h(v) =N(v)/D(v). (5) 


Letting N(v) be analytic except for a branch point 
at -{ and a pole at v,, and D(v) be analytic except 
for a branch point at the origin, one can immedi- 
ately arrive at the coupled integral equations 


xa ee dy Dw) D(v’) 


v'(v'- 





N(v) = (6) 


v)’ 


(v)2N(v") 
(v’-v,)(v’ - v) ° 





D(v)=1 ia Jaw (7) 


The subtraction is made in Eq. (7) to keep the 
consistent asymptotic behavior of h(v) — O[(1/v)Inv] 
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for f7#0, and h(v)~O(1/v) for f?=0. In the limit 
f?=0, the solution of Eqs. (6) and (7) reduces to 
the effective-range formula (1). The residue T 
is then trivially related to I’, of Eq. (3). Fora 
coupling constant greater than zero, Eqs. (6) and 
(7) can be solved by a straightforward iteration 
in (f?M). The series is uniformly convergent for 
(f?M) <1 (the actual radius of convergence may 
be greater than 1). In this calculation we use 

f? =0.08, (f?M) =0.533; I and y, are adjusted to 
fit two precisely known singlet S p-p phase shifts 
at 1.397 and 2.425 Mev,‘ or the deuteron binding 
energy and the triplet scattering length. Calcu- 
lated »“*cot5 curves are given in Figs. 1 and 2 
for the singlet p-p and the triplet n-p, respective- 
ly.5 Shape parameters have also been calculated, 
and the quadratic approximation 


vcoté =-1/a+4rv- Pri, (8) 


is plotted on the corresponding figures. 

We summarize our results as follows: 

(a) Since the effective “interaction pole” re- 
places only those branch points at »y <-1, the 
present result should be reasonably reliable up 
to v <1 (=40 Mev). 

(b) Although the power series in v diverges for 
|v| 2} because of the one-meson branch point, 
our calculated v“*coté curve remains quite close 
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FIG. 1. Singlet p-p scattering, with a=-5.5, r=2.0 
(44,=-2.6, P=0.073). 
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FIG. 2. Triplet n-p scattering, with a@=3.8, r=1.2 


(vy=-2.0, P=0.028). 


to a straight line up to y=1. 


(c) As f? goes from zero to 0.08, v, moves 
from -1.8 to -2.0 for the triplet amplitude and 
from -1.34 to -2.6 for the singlet p-p amplitude. 
The ratio (-I'/v,), which is an approximate meas- 
ure of the effect of the pole, decreases by 14% 
for the triplet and 30% for the singlet. These 
numerical results point to the fact that the one- 
meson force is far from sufficient to give the 
required attraction.* However, the qualitatively 
reasonable positions and strengths of these 
“interaction poles” have added to our confidence 
in replacing unknown branch cuts by such poles. 
Also, it indicates that the successive inclusion 
of outer branch cuts would make the role of the 
phenomenological “interaction pole” less and less 
important. A systematic procedure of successive 
approximations appears to be quite possible. 

We wish to thank Professor Geoffrey F. Chew 
for initiating this investigation and for many 


enlightening discussions. 





‘Work performed under the auspices of the U. S. 


Atomic Energy Commission. 


1§. Mandelstam, Phys. Rev. 112, 1344 (1958). Man- 
delstam has proved his representation through the sixth 
order in perturbation theory. J. Bowcock and D. 
Walecka (private communication) have proved the fixed- 
angle dispersion relation for the Born series of non- 
relativistic potential scattering with Yukawa-type poten 
tials. Independently, Blankenbecler, Goldberger, Khur', 
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and Trieman have proved the Mandelstam representa- 
tion for the same potential scattering problem to all 
orders in the Born series as well as the Fredholm 
series (private communication). 

*D. Y. Wong, Phys. Rev. Letters 2, 406 (1959). 
‘This method was suggested by G. F. Chew and S. 
Mandelstam, Lawrence Radiation Laboratory Report 
UCRL-8728, April 1, 1959 (unpublished), in connection 

with an analogous problem for pion-pion interaction. 
‘MacGregor, Moravesik, and Noyes, University of 
California Radiation Laboratory Report UCRL-5582-T; 


Knecht, Messelt, Berners, and Northcliffe, Phys. Rev. 


114, 550 (1959). 

“Sin making use of the two proton-proion phase shifts, 
we have simply assumed that v”*coté as calculated from 
Eqs. (6) and (7) is to be compared with the usual Cou- 
lomb modification 


{[2mqcoté/exp(2mm) - 1] + 2gnh(q)}; n=e?/fiv. 






SCini, Fubini, and Stanghellini; W. Alles and A. 
Tomasini; and S. Matsuyama (private communication) 
have made separate attempts to determine the coupling 
constant from dispersion relations using observed s- 
wave parameters. It is clear that at least three param- 
eters are needed for a reliable determination of the 
coupling constant. It seems to us that our present 
knowledge of the s wave is inadequate to give such a 
three-parameter set. However (Fubini and Stanghellini, 
private communication) if one assumes the value of 

f* known, the formulas for gcoté and its derivative as 
given by Cini, Fubini, and Stanghellini evaluated at 

¢ =-1/2 imply a positive shape parameter and a 'S, 
phase shift of ~ 48° at 40 Mev, as do our formulas. 
Preliminary results indicate that experiment may sup- 
port this conclusion [MacGregor, Moravcsik and Noyes, 
University of California Radiation Laboratory Report- 
UCRL-5582-T (unpublished)] . 





ERRATUM 





EVIDENCE FOR ANISOTROPY OF THE SUPER- 
CONDUCTING ENERGY GAP FROM ULTRA- 
SONIC ATTENUATION. R. W. Morse, T. Olsen, 
and J. D. Gavenda [Phys. Rev. Letters 3, 15 
(1959)]. 


Measurements reported in this paper for pro- 
pagation along the [100] direction have been 


found to be erroneously labelled. Actually these 
data were taken for a direction of propagation 
which was perpendicular to [001] and 18° from 
[100]. From subsequent measurements along 
[100] one estimates the limiting energy gap to be 
(3.540.1)kT, for this direction. 











VoLUME 3, NUMBER 4 


PHYSICAL REVIEW LETTERS 








AuGusT 15, 1959 














ABSTRACTS 


In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to be published)” rather than 
to this Journal. 





SIGNIFICANCE OF ELECTROMAGNETIC PO- 
TENTIALS IN THE QUANTUM THEORY. 

Y. Aharonov and D. Bohm, H. H. Wills Physics 
Laboratory, University of Bristol, Bristol, Eng- 
land (Received May 28, 1959; revised manuscript 
received June 16, 1959). 


In this paper, we discuss some interesting pro- 
perties of the electromagnetic potentials in the 
quantum domain. We show that, contrary to the 
conclusions of classical mechanics, there exist 
effects of potentials on charged particles, even 
in the region where all the fields (and therefore 
the forces on the particles) vanish. We then 
discuss possible experiments to test these con- 
clusions; and, finally, we suggest further pos- 
sible developments in the interpretation of the 
potentials. 


THEORY OF MULTIPLE SCATTERING: SECOND 
BORN APPROXIMATION AND CORRECTIONS 

TO MOLIERE’S WORK. B. P. Nigam, M. K. 
Sundaresan, and Ta-You Wu,* Division of Pure 
Physics, National Research Council, Ottawa, 
Canada (Received March 19, 1959). 


The formula given by Moliére for the scattering 
cross section of a charged particle by an atom, 
on which has been based the formula for the 
“screening angle” x a in his theory of multiple 
scattering, has been examined and found to 
contain an inconsistent approximation in all orders 
of the parameter a, =zZ/1378 except the lowest 
(the first Born approximation). In the present 
work, the correct expression of Dalitz is used 
for the single-scattering cross section of a 
relativistic Dirac particle by a screened atomic 
field up to the second Born approximation. It is 
found that the effect of the deviation from the 


first Born approximation on the screening angle 
is much smaller than Moliére’s expression for 
this quantity would lead one to believe. This is 
so because the deviation from the first Born 
approximation is very small at the small angles 
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—— 


that go into the definition of the screening angle, 
In Moliére’s work, all the effect of the deviation 
from the first Born approximation on the distri- 
bution function f(@) for multiple scattering is 
contained in the quantity B which depends only 
on x,- In the present work, it is shown that in 
a consistent treatment of terms of various 
orders in a,, there exist additional terms of 
orders zZ /137 in the distribution function. These 
terms, which represent the second Born approx- 
imation, become important at large angles. Cal- 
culations have been carried out for the scat- 
tering of 15.6-Mev electrons by Au and Be. The 
1/e widths of the distribution function obtained 
are in good agreement with the experimental 
result of Hanson et al., whereas Moliére’s 
theory gives too great a width compared with 
the experimental value for Be. 





- 
At present at the Institute for Advanced Study, 
Princeton, New Jersey. 


ELECTRIC FIELD DISTRIBUTIONS IN AN ION- 
IZED GAS. I. Michel Baranger and Bernard 

Mozer, Carnegie Institute of Technology, Pitts- 
burgh, Pennsylvania (Received March 12, 1959). 


A method for improving systematically the 
Holtsmark distribution is described. It is based 
on a cluster-type expansion and takes into ac- 
count increasing orders of correlation. In this 
paper, it is applied to the calculation of the dis- 
tribution of the high-frequency component of the 
electric field in an ionized gas in thermal equilib- 
rium. 


MOMENT OF INERTIA OF INTERACTING 
MANY-BODY FERMION SYSTEMS. R. D. Amado 
and K. A. Brueckner, University of Pennsylvania, 
Philadelphia, Pennsylvania (Received March 9, 
1959). 


It is shown that the moment of inertia of a non- 
interacting many-body fermion system moving 
under periodic boundary conditions has the clas- 
sical or rigid value when calculated on the 
“cranking” model of Inglis. By investigating the 
analogous “pushing” case for the inertial mass, 
we show that the rigid moment can be associated 
with rigid rotation in spite of apparent surface 
currents. The effect of particle-particle forces 
is investigated in the lowest order of perturba- 
tion theory. The terms corresponding to a level 
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shift or effective mass are just compensated by 
other terms and there is no change in the mo- 
ment of inertia. The possible general validity 
of these results and their consequence is dis- 
cussed. 


NEGATIVE TEMPERATURES. P. T. Landsberg, 
Department of Natural Philosophy, Marischal 
College, Aberdeen, Scotland (Received March 

6, 1959). 


It is shown by an example how the concept of 
negative temperatures finds a natural place in 
the approach to thermodynamics which was orig- 
inated by Carathéodory. 


COOLING BY ADIABATIC MAGNETIZATION. 
W.P. Wolf, Clarendon Laboratory, Oxford, Eng- 
land (Received March 20, 1959). 


A new method of producing temperatures below 
1°K is considered. It utilizes the fact that the 
entropy of certain paramagnetic salts is increased 
by the isothermal application of a magnetic field, 
so that conversely adiabatic magnetization should 
produce cooling. The factors limiting the temp- 
eratures which can be reached are discussed, and 
some suitable classes of salts proposed. 


MAGNETIC SUSCEPTIBILITY OF THORIUM 
METAL IN THE RANGE 130-300°K. J. F. Smith 
and J. D. Greiner, Institute for Atomic Research 
and Department of Chemistry, Iowa State Col- 
lege, Ames, Iowa (Received March 16, 1959). 


The magnetic susceptibility of thorium has 
been measured in the temperature range 130- 
$00°K. Within the precision of measurement, 
the susceptibility was found to remain constant 
at a value of + 0.410+0.002 emu/g. This sus- 
ceptibility behavior is in contrast to the elastic 
behavior which indicates a change in the char- 
acter of the interaction forces near 253°K. 


PARAMAGNETIC RESONANCE OF S-STATE 
IONS IN STRONTIUM CHLORIDE. W. Low and 

U. Rosenberger, Department of Physics, The 
Hebrew University, Jerusalem, Israel (Received 
February 20, 1959). 


Paramagnetic resonance spectra have been 
observed of Mn*+, Gd*+, and Eu+ in single crys- 





tals of strontium chloride grown from the melt. 
The Mn** shows a very small cubic field split- 
ting a <1x10™* cm™ and a hyperfine structure 
constant of A = 81.2x10~* cm™', indicating some 
covalent bonding. The Gd** and Eu** have cubic 
symmetry with the following splitting parameters: 
Gd**: c =+(39.6+0.1)x10™ cm™, d=#(0.2+0.1) 
x10~* cm“, g=1.9906+0.001; Eu?*: c =(52+5) 
x10 cm™, A! =(34.5+0.3) x10 cm™, AS 
=(15.5+0.3)x10™ cm™, g=1.995+0.001. The 
smaller cubic field splittings in SrCl, compared 
with CaF, seem to favor a mechanism responsible 
for the cubic field splitting which is linear in the 
crystal potential terms. 


ZEEMAN SPLITTINGS OF PARAMAGNETIC 
ATOMS IN CRYSTALLINE FIELDS. H. Statz, 
Raytheon Manufacturing Company, Waltham, 
Massachusetts, and G. F. Koster, Massachusetts 
Institute of Technology, Cambridge, Massachu- 
setts (Received February 24, 1959). 


The paper deals with the energy levels as a 
function of magnetic field of paramagnetic ions 
in crystalline surroundings. The necessary in- 
formation is given to write down the Hamiltonian 
matrix for any ion in any crystalline field. The 
method has been previously developed and it 
consists of determining the Hamiltonian matrix 
of 8H-(L+ 28) for the lowest states of an ion by 
the methods of group theory and time-inversion 
considerations. The Hamiltonian matrix contains 
a certain number of undetermined constants 
which may be evaluated if the wave functions of 
the ions are known, or the constants may be de- 
termined from experiment. In general, the num- 
ber of available constants is larger than found 
in the conventional spin Hamiltonian formalism. 
The present method may be applied in cases 
where the conventional spin Hamiltonian for- 
malism is inadequate of describing experimental 
results or where greater accuracy is required. 


NERNST AND ETTINGSHAUSEN EFFECTS IN 
GERMANIUM BETWEEN 300 AND 700°K. Herbert 
Mette, Wolfgang W. Gartner, and Claire Loscoe, 
U. S. Army Signal Research and Development 
Laboratory, Fort Monmouth, New Jersey (Re- 
ceived March 6, 1959). 


The Nernst and Ettingshausen effects in ger- 
manium single crystals of different conductivity 
type and with various impurity densities have 
been measured between 300 and 700°K in magnetic 
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fields of 9000 gauss, the Nernst effect also at 
2100 gauss. The experimental results are com- 
pared with theoretical expressions for the Nernst 
and Ettingshausen coefficients and with previously 
reported values for the thermal conductivity, 
which links the two effects through the Bridgman 
relationship. The qualitative agreement is very 
good; quantitative discrepancies are explained in 
terms of deviations of the sample mobilities 
from the theoretically assumed pure low-field 
lattice mobilities and by uncertainties in the 
numerical values used for the theoretical calcu- 
lations. A derivation of the Nernst constant in 
the range of negligible phonon-drag effect is 
given, based on the analogy between thermo- 

and photomagnetoelectric effects. 


DISLOCATION DYNAMICS AT LOW TEMPERA- 
TURES. J. Lothe, Universitetets Fysiske Insti- 
tutt, Blindern, Oslo, Norway, and J. P. Hirth, 
Metals Research Laboratory, Carnegie Institute 
of Technology, Pittsburgh, Pennsylvania (Re- 
ceived March 12, 1959). 


A model for low-temperature dislocation motion 
is developed using the concept of nucleation and 
growth of kinks in dislocations lying in Peierls’ 
potential troughs. The model is compared with 
existing experimental data. 


STATISTICAL MECHANICAL THEORY OF A 
RANDOM FERROMAGNETIC SYSTEM. R. Brout, 
Department of Physics, Cornell University, 
Ithaca, New York, and Bell Telephone Labora- 
tories, Murray Hill, New Jersey (Received 
March 11, 1959). 


The behavior of solid solutions of paramagnetic 
impurities which are exchange-coupled in a non- 
magnetic substrate turns out to yield a consider- 
able body of information with regard to the nature 
of the exchange coupling as well as detailed tem- 
perature dependence of the spin system. In the 
present paper, a rigorous expansion of the mean 
free energy averaged over random sites is pre- 
sented. It is shown that a ferromagnetic phase 
transition does occur. The Curie point is given 
as a function of concentration for the case of 
weak dilution in an implicit power series form. 

Many interesting qualitative features arise in 
the study of these systems. If the curve of mag- 
netic moment vs temperature has inflections 
this indicates short-range exchange forces, 
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whereas smooth curves indicate long-range for- 
ces. Similarly, long-range forces give rise to 
smooth behavior of the Curie point as a function 
of concentration for dilute samples. Alternatively, 
short-range forces give rather violent changes 
in Curie point near atomic fractions ~ (number of 
nearest neighbors)~*. The method of series de- 
velopment used in this paper gives rise to this 
expected qualitative behavior and also enables 
one to make quantitative prediction if the ex- 
change potential is known. 

Consideration is also given to the antiferro- 
magnetic analog together with a discussion of 
expected behavior of such systems in a resonance 
experiment. 


APPLICATION OF THE DIFFUSION -MODIFIED 
BLOCH EQUATION TO ELECTRON SPIN RESO- 
NANCE IN ORDINARY AND FERROMAGNETIC 
METALS. Jerome I. Kaplan, United States Naval 
Research Laboratory, Washington, D. C. (Re- 
ceived March 12, 1959). 


A new solution to the electron spin resonance 
absorption in metals is given by using Bloch’s 
equation, modified to include diffusion. It is 
shown that the coefficient for a similar diffusion 
term added to the Landau-Lifshitz equation for 
ferromagnetic resonance must be of order 10* 
times smaller than that for nonmagnetic metals. 
A restricted formula for saturation in nonmagnetic 
metals is given. 


OPTICAL PROPERTIES OF ACTIVATED AND 
UNACTIVATED HEXAGONAL ZnS SINGLE 
CRYSTALS. S. P. Keller and G. D. Pettit, Re- 
search Laboratory, International Business Ma- 
chines Corporation, Poughkeepsie, New York 
(Received March 13, 1959). 


Single crystals of hexagonal ZnS, unactivated 
and activated with Cu, Al, or Mn, have been ex- 
amined. The polarization effects in the trans- 
mission spectrum of unactivated ZnS were meas- 
ured further into the fundamental absorption 
than were previously measured. The theoretical 
prediction that, for direct transitions at k=0, 
light polarized perpendicular to the c axis is 
more strongly absorbed than light polarized 
parallel to the c axis is borne out throughout the 
fundamental absorption region except in the 
wavelength region between 290 and 325 mu. The 
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wavelength dependence and the polarization of 
the excitation and fluorescence spectra of the 
activated and unactivated crystals were meas- 
ured at room temperature and at 77°K. The ex- 
citation spectra showed an agreement with the 
selection rule at the edge of the fundamental 
absorption, but there was a reversal of the se- 
lection rule deep in the absorption region. 

There was also an impurity absorption exhibiting 
the same polarization properties as the edge. 
Some fluorescent emissions were polarized per- 
pendicular and some were polarized parallel to 
the c axis. Speculations are made on the re- 
versal of the polarization deep in the fundamental 
absorption and on the symmetry and the nature 

of the sites causing the various fluorescence 
bands. 


POLARIZATION PROCESSES IN PARTIALLY 
ILLUMINATED PHOTOCONDUCTING INSULA- 
TORS. M. Y. Ben Sira and B. Pratt, Technion, 
Israel Institute of Technology, Haifa, Israel, 
and E. Harnik and A. Many, Department of 
Physics, Hebrew University, Jerusalem, Israel 
(Received March 9, 1959). 


A quantitative analysis of polarization phenom - 
ena in photoconducting insulators is presented 
for the case in which a strip of the sample is 
illuminated far from the electrodes. It is shown 
that if the drift length of the photoexcited car- 
riers is small compared with the dimensions of 
the dark and illuminated regions, the photocur- 
rent decays exponentially during both the buildup 
and release of polarization. The decay constant 
is essentially independent of the characteristics 
of the contacts and the applied voltage, and can 
be expressed in terms of geometrical factors 
and the photoconductance of the illuminated strip. 
The conclusions of the analysis are checked in 
colored single crystals of NaCl. The experi- 
mental results are found to be in satisfactory 
agreement with the theory presented. 


KINETICS OF VACANCY MOTION IN HIGH- 
PURITY ALUMINUM. Warren DeSorbo and David 
Turnbull, General Electric Research Laboratory, 


Schenectady, New York (Received March 12, 
1959). 


It is shown that in very-high-purity aluminum, 
Aly the rate of vacancy annealing depends on 








vacancy concentration and annealing temperature 
but is independent of the temperature T; of va- 
cancy injection per se. The rate can be described 
as the sum of first and second order components. 
The first order component becomes most pro- 
minent at a monovacancy concentration estimated. 
to be 10~® atom fraction. It is shown that the re- 
sults are consistent with the Koehler-Seitz- 
Bauerle dissociative mechanism. The activation 
energy for diffusion of monovacancies in Aly is 
found to be 0.65+0.06 ev. This, combined with 
earlier results on the formation energy of vacan- 
cies, gives 1.44+0.11 ev for the activation en- 
ergy for self-diffusion in aluminum by a mono- 
vacancy mechanism. 

In zone-refined aluminum, Al,, of lesser pur- 
ity, the rate of vacancy annealing depends upon 
T; per se and falls off more rapidly with de- 
creasing vacancy concentration than in Aly: Two 
hypotheses for the impurity effects are consid- 
ered, namely: (1) trapping of vacancies by im- 
purity atoms and (2) inhibition of dislocation 
climb by adsorbed impurities. 


EFFECT OF OXYGEN ON THE WORK FUNCTION 
OF BARIUM. P. A. Anderson and A. L. Hunt, 
State College of Washington, Pullman, Washing- 
ton (Received March 2, 1959). 


While a high sensitivity of the work function to 
gaseous contaminants is generally characteristic 
of metals which adsorb the common gases 
strongly, barium combines outstanding work 
function stability with an exceptional affinity for 
these gases. The paradox is investigated by 
administering measured quantities of oxygen to 
vapor -deposited barium surfaces at known rates, 
with concurrent monitoring of the work function. 
It is found (1) that when clean barium is exposed 
suddenly to a massive dose of oxygen its work 
function is lowered by 0.31 ev; (2) that the work 
function of the heavily oxygenated surface re- 
mains constant for 5 days and then drifts very 
slowly toward the clean-metal value, and (3) 
that when oxygen is administered to a barium 
surface slowly (2x10"* molecules per cm? per 
hour) no change of work function occurs during 
or after the absorption of a quantity of oxygen 
equivalent to 100 complete monolayers of the 
oxide. 

The results are interpreted as stemming from 
two processes which compete for control of the 
surface: the expected surface oxidation and a 
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restorative process involving engulfment of 
oxygen ions. The mechanism of engulfment is 
discussed. 


WORK FUNCTION OF GOLD. Paul A. Anderson, 
State College of Washington, Pullman, Washing- 
ton (Received March 2, 1959). 


The work functions of sixteen gold surfaces 
have been determined by measurement of their 
contact differences of potential with respect to 
barium reference surfaces of known work func- 
tion. The gold surfaces were prepared by sub- 
jecting spectroscopic standard gold to repetitive 
fusion-solidification outgassing followed by frac- 
tional distillation and condensation on tantalum; 
the barium surfaces by a similar technique 
which yields surfaces reproducible and constant 
to 0.01 ev. Measurement was by the Kelvin 
method with a time interval of 15 seconds bet- 
ween deposition and measurement of a fresh 
surface. 

The contact difference of potential Au-Ba is 
found to be 2.31+40.02 volts and the work function 
of gold 4.83+0.02 ev, referred to a barium work 
function of 2.52 ev. 


THEORY OF THE NUCLEAR MAGNETIC RESO- 
NANCE SHIFT IN PARAMAGNETIC CRYSTALS. 
F. Keffer, T. Oguchi,* W. O’Sullivan,f and J. 
Yamashita,* University of Pittsburgh, Pittsburgh, 
Pennsylvania (Received October 29, 1958; re- 
vised manuscript received April 27, 1959). 


A theoretical study is made of the shift of the 
F’* nuclear resonance in paramagnetic and anti- 
ferromagnetic MnF, which has been observed by 
Shulman and Jaccarino. The problem is reduced 
to that of a single Mn-F pair. A net hyperfine in- 
teraction is shown to arise from overlap effects in 
the ground-state ionic configuration (3d)5(2s)?(2p)® 
and from overlap and transfer effects to the con- 
figurations (3d)*°(2s)(2p)* and (3d)°(2s)?(2p)®°. These 
three configurations are equivalent to a single 
configuration involving bonding-type molecular 
orbitals. The results are in reasonable agree- 
ment with the experiment, the theoretical iso- 
tropic shift being slightly too small and the theo- 
retical anisotropic shift (small non-dipolar part) 
being slightly too large. A reappraisal is made 
of Tinkham’s data on paramagnetic resonance of 
Mn** impurities in ZnF,, which Bleaney has 
shown to be closely related to the Shulman- 
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Jaccarino data. It is found that there is no need 
to include, as did Tinkham, a large fraction of 
fluorine 3s and 3p functions into the bond. 


“On leave from Tokyo University of Education, 
Tokyo, Japan. 
t Present address: U. S. Naval Postgraduate School, 
Monterey, California. 
Present address: Institute for Solid State Physics, 
University of Tokyo, Tokyo, Japan. 


MEISSNER-OCHSENFELD EFFECT IN THE 
BOGOLJUBOV THEORY. R. M. May and M. R. 
Schafroth,f The F. B. S. Falkiner Nuclear Re- 
search and Adolph Basser Computing Labora- 
tories, School of Physics, The University of 
Sydney, Sydney, N.S.W., Australia (Received 
January 19, 1959). 


The magnetic response of the Bogoljubov 
model is calculated, by employing a new tech- 
nique for dealing with gauge-noncovariant ap- 
proximations to gauge-covariant Hamiltonians. 

By extending our treatment to very high orders 
of perturbation theory, we derive a Meissner- 
Ochsenfeld effect with infinite correlation length, 
as anticipated by London. For larger wave- 
lengths, the response differs from the London 
one, in the way predicted by Pippard. 


"peceased " 


CRYSTAL ATOM SPACING AND THE WAVE 
LENGTH OF AN INCIDENT ATOM BEAM. 
Robert E. Davidson, Langley Aeronautical Lab- 
oratory, National Advisory Committee for Aero- 
nautics, Langley Field, Virginia (Received 
August 1, 1958; revised manuscript received 
February 6, 1959). 


A beam of atoms is assumed to be incident on 
a crystal, and a certain aspect of the energy 
exchange is investigated insofar as is possible 
without elaborate experiments or calculations. 
The interaction potential is assumed to be that 
for helium or lithium fluoride and is known from 
other sources. By taking account of the thermal 
motion of the crystal atoms tangentially to the 
surface of the crystal, it is deduced that the 
probability that the beam will give up energy to 
the surface becomes relatively large when the 
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de Broglie wavelength of the tangential component 
of beam momentum is decreased to the same 
order of magnitude as the spacing of atoms in the 
crystal. 


CONVERSION ELECTRON ANGULAR CORRELA- 
TIONS: GENERAL K-SHELL FORMULATION 
AND THRESHOLD LIMIT. Robert C. Young, The 
Rice Institute, Houston, Texas (Received March 
13, 1959). 


A general treatment of polarization-angular 
correlations involving conversion electrons is 
presented. Formulas are given for the b,, coef- 
ficients for all types of polarization of the con- 
version electron, electric and magnetic multi- 
poles, and arbitrary complexity v. It is shown 
that only three coefficients are needed to com- 
pletely characterize correlations of conversion 
electrons from the K shell, for a given type of 
transition. Curves are given of the threshold 
values for these coefficients for L=1-5, electric 
and magnetic, versus atomic number. The re- 
lativistic Coulomb functions (point nucleus, no 
screening) are presented, along with a number 
of properties and expansions, in a convenient 
and consistent notation. Several identities and 
relations involving the vector addition coefficients 
are included. 


MOST PROBABLE ENERGY LOSS OF FAST 
ELECTRONS. H. E. Hall, A. O. Hanson, and 
D, Jamnik, University of Illinois, Urbana, Illi- 
nis (Received March 6, 1959). 


The most probable energy losses of electrons 
of 3.4, 7.2, 11.2, 15.6, and 19.4 Mev have been 
measured for Al, Cu, Ag, Ta, Au, and U sam- 
ples of about 1 g/cm’. The measurements agree 
with available calculations and refute the previ- 
ous discrepancy for gold. 


NEUTRON POLARIZATION IN THE REACTION 


D(d,n)He* AT E,=8.2 Mev. Wilfried W. Daehnick,” 


Wayman Crow Laboratory of Physics, Washing- 
ton University, St. Louis, Missouri (Received 
March 16, 1959). 


The polarization of neutrons from the reaction 
D(dn)He* was measured at 8.2 Mev incident 
deuteron energy. By the use of He* as a polariza- 
tion analyzer and a triple coincidence technique 


systematical errors in the experiment were kept 
very small, but at a sacrifice in counting rate. 
Polarization measurements were made at center- 
of-mass angles of 0°, 47°, and 59°. The values 
obtained for the neutron polarization and their 
probable statistical errors are P(0) =(-)0.007 
+0.007, P(47)=-0.097+0.060, and P(59) =-0.101 
+0.060. 


a 
Now at Palmer Physical Laboratory, Princeton 
University, Princeton, New Jersey. 


EVIDENCE FOR THE POLARIZATION OF B” 
NUCLEI PRODUCED IN A (d,p) REACTION. 

L. F. Chase, Jr.,* and G. Igo,? Stanford Univer- 
sity, Stanford, California (Received March 2, 
1959). 


Evidence is presented that a lower limit for 
the polarization of B’* nuclei produced in the 
(d,p) reaction at 2.8 Mev is 17+5%. The sign 
of the polarization Pp is in disagreement with 
the Newns model, and it is in agreement with a 
modified classical model. The sign of the pro- 
ton polarization determined by Hillman and by 
Juveland and Jentschke is also in agreement 
with the modified classical model. 


*Now at Lockheed Aircraft Corporation, Missile Sys- 
tems Division, Palo Alto, California. 

TNow at the University of Heidelberg, Heidelberg, 
Germany. 


(p,n) CROSS SECTIONS AND PROTON OPTICAL- 
MODEL PARAMETERS IN THE 4- TO 5.5-Mev 
ENERGY REGION. Richard D. Albert, Univer- 
sity of California, Lawrence Radiation Labora- 
tory, Livermore, California (Received March 

23, 1959). 


(p,m) cross sections were measured for 18 
medium-weight nuclei in the 4- to 5.5-Mev re- 
gion. In this energy region, de-excitation of the 
compound nucleus by neutron emission is largely 
favored over charged-particle emission; con- 
sequently, the (p,m) reaction dominates com- 
peting compound-elastic and inelastic scattering 
processes. The experimental results are com- 
pared with theoretical reaction cross sections 
predicted by the optical model for protons. It is 
found that the parameters which provide the best 
theoretical fit to these data also provide reason- 
able fits to neutron total, nonelastic, and elastic 
cross sections measured elsewhere when the 
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latter are compared with theoretical cross sec- 
tions predicted by the optical model for neutrons. 
It is concluded that the intrinsic nuclear well 
depths for protons and neutrons are equal to 
+7% uncertainty. 


DECAY OF Pm™. Kenneth S. Toth and O. B. 
Nielsen, Institute for Theoretical Physics, Uni- 
versity of Copenhagen, Copenhagen, Denmark 
(Received March 6, 1959). 


Pm‘ has been produced by an (a,) reaction 
on Pr**! in the Copenhagen cyclotron. The half- 
life of Pm™ has been determined and found to be 
closer to 450 days rather than to the previously 
reported value of 300 days. The decay of the 
nuclide to levels in Nd’ has been studied, using 
a 100-channel gamma analyzer and a six-gap 
8-ray spectrometer. Three y rays have been 
found, with the following energies: 474, 615, 
and 695 kev, respectively. It has been determined 
that the three y rays are E2 transitions and that 
they are in cascade with one another. A decay 
scheme for Pm™, consistent with the experi- 
mental data, has been proposed. 


708-kev GAMMA-RAY IN P*®. J. J. Singh, De- 
partment of Physics, University of Kansas, 
Lawrence, Kansas (Received December 4, 1958; 
revised manuscript received May 14, 1959). 


The purpose of this report is to present evidence 
that the 708-kev gamma ray in P®, arising from 
a transition between the second excited state and 
the ground state, is a strongly mixed M1 and E2 
transition. It is suggested that the observed 
E2/M1 amplitude ratio indicates the operation 
of selection rules for M1 radiation in self-con- 
jugate nuclei. 


SPONTANEOUS NUCLEAR REACTION IN MOLEC- 
ULAR HYDROGEN. Richard G. Brewer,* Depart- 
ment of Chemistry, University of California, 
Berkeley, California (Received March 4, 1959). 


The probability that two protons in a hydrogen 
molecule spontaneously react to form a deuteron 
has been calculated by the use of WKB wave 
functions. The complete shape of the molecular 
barrier has been derived from spectroscopic 
data for the hydrogen molecule and from a per- 
turbation calculation of the helium atom. Use of 
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this potential function showed that the probability 


of penetrating the molecular barrier (for a groun- 


state molecule) is greater by a factor of 10" than 
of penetrating a Coulomb barrier at the same 
internuclear distance. The half-lives for nuclear 
reaction in the isotopic molecules H*H’, H*H?, 
and H?H? are 5x10, 5x10**, and 1x10™ years, 


* 
Present Address: Department of Chemistry, 
Harvard University ,Cambridge, Massachusetts. 


DETERMINATION OF THE DIPOLE MOMENT 
AND ISOTOPE SHIFT OF RADIOACTIVE Hg’ 
BY “DOUBLE RESONANCE.” Adrian C. Melissi- 
nos,* Department of Physics and Research La- 
boratory of Electronics, Massachusetts Institute 
of Technology, Cambridge, Massachusetts (Re- 
ceived January 8, 1959). 


Paramagnetic resonance was established be- 
tween m -sublevels (Am =+1, AF =0) of the *P, 
state of radioactive Hg’®” at 3000 Mc/sec. From 
these data the nuclear interaction constant A 
was found to be (513.541) x10-* cm™?, and bar- 
ring hfs anomalies it led to a ratio of moments 
Lyor/Hy99 =A y97/A joo = 1-045; further, the nuclear 
spin of Hg’” was ascertained to be 3. The 
double resonance was combined with magneto- 
optic scanning to give the isotope shift of Hg’, 
which was found to be in the 2537 A line +(91+5) 
x 107° cm~! from Hg’™®. The radioactive mer- 
cury was produced by the Au’®"(d, 2m)Hg'™ reac- 
tion and used in vapor form. Satisfactory sig- 
nals were obtained with as few as 3x10” atoms. 


“present address: Department of Physics, Univer- 
sity of Rochester, Rochester, New York. 


DIPOLE AND QUADRUPOLE MOMENTS OF THE 
ISOMERIC Hg’*’** NUCLEUS; ISOMERIC ISOTOPE 
SHIFT. Adrian C. Melissinos,* Department of 
Physics and Research Laboratory of Electronics, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts, and Sumner P. Davis, 
Department of Physics and Spectroscopy Labora- 
tory, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received January 8, 
1959). 


The hyperfine structures of five optical lines of 
65-hour radioactive Hg” and 25-hour isomeric 
Hg’*’* were measured. Radioactive mercury pro 
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electrodeless discharge tubes containing as few 
as 5X10*” atoms, free of natural mercury, gave 
adequate light intensities. The analysis of the 
hyperfine structure confirmed the magnetic di- 
pole moment and isotope shift of Hg’®” in its 
ground state (I=4), and gave the following values 
for Hg’”* (based on J=13/2): pyg7x = -1.0420.01 
nm; Qiore = (1.540.3) x10°™™* cm?; isotope shift 
displacement (from Hg’ in the 2537 A line) 
(,70+7)x10-* cm™*. Thus an isomeric shift of 
atomic energy levels resulting from the excita- 
tion of the nucleus from its ground state J=4 to 
the isomeric state J= 13/2 is observed. In the 
2537 A line of mercury this shift amounts to 
(-21+6) x10-* cm™. 

* Present address: Department of Physics, Univer- 
sity of Rochester, Rochester, New York. 


RADIOCHEMICAL STUDIES OF THE FISSION 
OF U** INDUCED BY HELIUM IONS. R. Gunn- 
ink and J. W. Cobble, Department of Chemistry, 
Purdue University, Lafayette, Indiana (Received 
March 9, 1959). 


The absolute fission yields of approximately 
twenty-five nuclides from the helium-ion- induced 
fission of U*** were determined with an accuracy 
of 5-15% for several energies ranging from 20 
to40 Mev. Such features of fission as the sym- 
metric-asymmetric modes of fission, the rela- 
tion of total fission cross section to compound 
nucleus theory, fine structure in fission product 
distribution, valley to peak ratios, and neutron 
emission are discussed as well as some of the 
experimental detail involved. 


SEMICLASSICAL THEORY OF INTERNAL RAY- 
LEIGH SCATTERING. A. M. Cormack,* Physics 
Department, University of Cape Town, Ronde- 
bosch, Cape, South Africa (Received May 15, 
1958; revised manuscript received April 6, 

1959), 


The angular distribution of the gamma rays 
emitted by a bare nucleus with a fixed orienta- 
tion in space may be modified if the nucleus is 
surrounded by electrons, because of the elastic 
scattering of the gamma rays by the electrons. 
Ageneral expression is given for the angular 
distribution of the radiation from arbitrary elec- 
tric and magnetic multipoles scattered by ar- 
bitrary electron distributions on the basis of a 
classical model, and an “internal scattering 


coefficient” is given explicitly for spherically 
symmetrical electron distributions. Although 
the radiation scattered from spherically symmet- 
rical charge distributions is indistinguishable 
from the unscattered radiation, the internal scat- 
tering coefficient has been calculated for these 

on the basis of a Thomas-Fermi atom so that an 
indication may be obtained of the smaliness of 
the scattering from asymmetries in atomic elec- 
tron distributions due to molecular or crystal 
binding effects. 


“Now at Tufts University, Medford, Massachusetts. 


ELECTRON PAIR PRODUCTION IN THE FIELD 
OF THE PROTON AND IN THE FIELD OF THE 
ELECTRON BY PHOTONS OF ENERGY FROM 
10 Mev TO 1 Bev. E. L. Hart, G. Cocconi, V. T. 
Cocconi, and J. M. Sellen,* Cornell University, 
Ithaca, New York (Received March 12, 1959). 


A 24-in. diffusion cloud chamber filled with 
hydrogen and located in a magnetic field has been 
placed in the hardened bremsstrahlung beam of 
the Cornell synchrotron to study electron pair 
production in the proton field (PF-pairs) and in 
the electron field (EF-pairs). The EF-pairs 
could be detected with an average efficiency of 
~ 85 % and are about as abundant as the PF-pairs. 
A total of 3065 pairs produced by photons with 
energy from 10 to 1040 Mev have been analyzed. 

The principal results are: (1) The recoil mo- 
mentum distribution of the EF-pairs is in good 
agreement with the distribution predicted by Suh 
and Bethe. (2) The cross section for EF-pairs 
is consistent with that calculated by Wheeler and 
Lamb. (3) The distribution of the opening angle 
of PF- and EF -pairs is in general agreement 
with Borsellino’s calculations, at large angles. 
However, the peak appears at smaller angles 
than predicted. 


o 
Now at the Ramo-Wooldridge Corporation, Los 
Angeles, California. 


INTERSECTING BEAM SYSTEMS WITH STOR- 
AGE RINGS. G. K. O’Neill and E. J. Woods, 
Palmer Physical Laboratory, Princeton Univer- 
sity, Princeton, New Jersey (Received Novem- 
ber 24, 1958). 


The equivalence of fixed- and variable-field 
particle acceleration systems for the adiabatic 
damping of synchrotron oscillations is pointed 
out. These two quite different acceleration 
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methods are therefore able to produce particle 
beams of equal density for beam stacking pur - 
poses. 

The transfer mechanism between an acceler- 
ator and a storage ring is discussed, and the 
properties of a fast-rise 3-kilogauss beam- 
switching magnet are shown. It is concluded 
that the source interaction rate obtainable in a 
proton storage ring system would be nearly 
independent of the focusing properties and repe- 
tition rate of the injecting accelerator. 

An improved design for intersecting-beam 
storage rings is described, in which several 
well-separated interaction regions could be used 
for simultaneous experiments. Standard types of 
alternating-gradient magnets would be required, 
and the over-all weight of synchrotron plus stor- 
age rings would be about one tenth as large as 
that of comparable beam-stacking accelerators. 


RECOIL MOMENTUM DISTRIBUTION IN ELEC- 
TRON PAIR PRODUCTION. K. S. Suh and H. A. 
Bethe, Laboratory of Nuclear Studies, Cornell 
University, Ithaca, New York (Received March 
12, 1959). 


Electron pair production by a very energetic 
photon in the field of a particle of arbitrary mass 
(in particular in the fields of an electron and a 
nucleus) is studied following the work of Borsel- 
lino. The distribution of recoil momenta gq is 
calculated for g of order of the electron mass, 
and it is shown that the recoil distribution is in- 
dependent of the mass of the recoil particle if 
appropriate variables are used. It is also ex- 
plicitly shown that the mass of the recoil particle 
does not make any difference in the recoil distri- 
bution for very small g (of order q,,;,)- The 
total cross section must therefore be independent 
of the mass of the recoil particle in the high- 
energy limit, as previously stated by Borsellino. 
The Wheeler-Lamb result for pair production in 
the field of a bound electron is also justified. 
The results also describe the electromagnetic 
production of any fermion pair if certain restric- 
tions are satisfied. 


ABSORPTION OF 1-Bev PHOTONS. E. Malamud, 


Laboratory of Nuclear Studies, Cornell Univer- 
sity, Ithaca, New York (Received March 12, 
1959). 


The total cross section for the attenuation of 
high-energy gamma rays has been measured in 
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various elements. The variation with atomic 
number, Z, has been observed by measuring 
the absorption of 1-Bev gamma rays in 12 dif- 
ferent elements ranging from hydrogen to ura- 
nium. Additional measurements were made in 
copper at 400 and 700 Mev to show the energy 
dependence of the absorption processes. These 
results are then combined with those of other 
investigators at lower energies and show that 
the theory of pair production in the nuclear field 
correctly predicts the cross section as a function 
of atomic number and photon energy. The meas- 
urements in low-Z elements give information on 
pair production in the field of the electron. The 
results are in closer agreement to the calcula- 
tions of Wheeler and Lamb than to the estimate 
of Joseph and Rohrlich. In addition, a short ex- 
periment was performed to measure the symme- 
try of the energy distribution between the elec- 
tron and positron members of the pair. The re- 
sults show that the energy-sharing curve is 
symmetrical as predicted by theory. 


PARITIES OF K-HYPERON PAIRS. Richard H. 
Capps, Laboratory of Nuclear Studies, Cornell 
University, Ithaca, New York (Received March 
13, 1959). 


The possibility of determining the relative 
parity of K-hyperon pairs by measurements of 
reactions of the type 7(or K)+p—K(or 1) + hyperon 
is investigated. If the protons are polarized ina 
direction perpendicular to the incident meson 
beam, it is shown that measurements of the cor- 
relation between the proton spin direction and 
the direction of the meson from the hyperon de- 
cay are capable of determining the parity un- 
ambiguously, provided that the K spin is zero 
and the hyperon spin is known. Several alter- 
nate parity-determining experiments that do not 
require polarized proton targets are discussed; 
these experiments are all quite difficult, how- 
ever. 


DETECTION AND GENERATION OF GRAVITA- 
TIONAL WAVES. J. Weber, University of Mary- 
land, College Park, Maryland (Received Feb- 
ruary 9, 1959; revised manuscript received 
July 20, 1959). 


Methods are proposed for measurement of the 
Riemann tensor and detection of gravitational 
waves. These make use of the fact that relative 
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motion of mass points, or strains in a crystal, 
can be produced by the space derivatives of the 
gravitational fields. The strains in a crystal 

may result in electric polarization in consequence 
of the piezoelectric effect. Measurement of volt- 
ages then enables certain components of the 
Riemann tensor to be determined. Mathematical 
analysis of the limitations is given. Arrange- 
ments are presented for search for gravitational 
radiation. 

The generation of gravitational waves in the 
laboratory is discussed. New methods are pro- 
posed which employ electrically induced stresses 
in crystals. These give approximately a seven- 
teen order increase in radiation over a spinning 
rod of the same length as the crystal. At the 
same frequency the crystal gives radiation which 
is about thirty-nine orders greater than that of 
a spinning rod. 


PROTON-NEUTRON MASS DIFFERENCE. 

Sigenobu Sunakawa* and Katsumi Tanaka, Argonne 
National Laboratory, Lemont, Illinois (Received 
March 6, 1959). 


An expression for the electromagnetic self- 
mass of a physical nucleon is obtained from the 
basic postulates of field theory. A complete set 
of intermediate states is then introduced and the 
one-nucleon state and the state with one nucleon 
plus a nucleon pair are taken into account. Direct 
comparison is thus established between the vari- 
ous methods of obtaining the proton-neutron mass 
difference. The mass difference is computed by 
two methods which use the Stanford form factors. 
It is shown that if the form factor which agrees 
in the experimental region is used in the non- 
experimental region in the nucleon-pair term, 
then there is a disagreement with the present 
electromagnetic explanation of the proton-neutron 
mass difference. 


* 
On ieave from Osaka University, Osaka, Japan. 


TIME RETARDATION IN STATIC AND STATION- 
ARY SPHERICAL AND ELLIPTIC SPACES. John 
Kronsbein and Erich A. Farber, University of 
a Gainesville, Florida (Received March 

, 1959). 


The well-known clock paradox is shown to be 
resolvable and consistent within the framework 
of the Lorentz group, using fewer assumptions 


about the observers than previously. In particu- 
lar it is not necessary in the spaces indicated to 
consider the effects of initial acceleration and 
final deceleration of the observers’ clocks since 
the existence of certain parallel periodic orbits 
makes it possible to confine observations to times 
of arbitrary length during which the relative ve- 
locity remains constant in a sense previously 
described. Complete agreement between the 
observers is attained, thus showing that despite 
the notion of “relative motion” it is possible to 
decide which of the two observers is subject to 
time retardation. 


NONLINEAR EFFECTS OF GRAVITATIONAL 
RADIATION. Asher Peres and Nathan Rosen, 
Department of Physics, Israel Institute of Tech- 
nology, Haifa, Israel (Received January 22, 
1959; revised manuscript received May 21, 
1959). 


It is shown that the cumulative effects of the 
nonlinear terms in Einstein’s equations rule out 
the possibility of stable small oscillations of a 
gravitational field about an equilibrium state, if 
the latter is supposed to be Minkowskian at in- 
finity. The perturbation, if unlimited in time, 
rather tends to take infinite values at large dis- 
tances from its source. This can be interpreted 
as an instability of gravitational radiation fields, 
and raises doubts concerning the validity of the 
linearized theory at large distances from the 
source. 


DISPERSION RELATION FOR NONLINEAR VAC- 
UUM POLARIZATION EFFECTS. J. G. Gilson, 
Chelsea College of Science and Technology, 
Manresa Road, London, England (Received 
March 9, 1959). 


The current induced in the vacuum by an ex- 
ternal field is a nonlinear functional of the latter 
when one goes beyond the first approximation. 

It is demonstrated that when this induced cur- 
rent is suitably defined it will satisfy the cau- 
sality requirement of being zero before the ex- 
ternal field is switched on. From this causality 
principle an exact dispersion relation is deduced. 
This relaxation of the linearity restriction in 
deriving a dispersion relation from causality 
implies that it may be possible to drop the lin- 
earity restriction in discussing the equivalence 
of causality and dispersion relations for other 
systems. 
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